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Abstract 

The  objective  of  this  report  is  to  identify  representations  and  issues  for  the  generic  assembly 
information  model  and  its  use  in  a proactive  tolerance  synthesis  scheme.  It  is  proposed  to  use  the  small 
displacement  torsors  (the  screw  parameters)  as  a mathematical  tool  for  analyzing  the  causal 
relationship  between  variations  in  parts’  shapes  and  sizes  (due  to  the  small  allowances  permitted  by 
tolerance  specifications),  and  their  effects  on  relative  positioning  on  adjacent  parts.  For  a given 
assembly,  the  allowable  geometric  variations  of  assembled  parts,  the  assembly  specification  and  the 
desired  functional  requirements  have  been  characterized  by  finite  sets  of  equality  and  inequality 
constraints.  An  optimization  process  is  then  formulated  to  study  the  tolerance  synthesis  and  analysis 
problems  in  the  assembly. 
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1.  Introduction 

This  report  gives  details  of  work  on  tolerance  synthesis  subsequent  to  our  earlier  work  [ROYOl],  For 
tolerance  synthesis  and  analysis,  we  principally  need  a detailed  description  of  the  “kinematic 
functions”  of  the  assembly,  by  which  we  mean  those  functions  defined  essentially  by  the  location,  size 
and  shape  (form)  of  associated  mating  features.  These  kinematic  functional  specifications  are  not 
directly  provided  by  the  customer  or  by  early  specifications  of  the  desired  product  / assembly  function. 
They  are  slowly  evolved  with  the  assembly  as  the  design  takes  concrete  shape  and  size  in  the  later 
phases  of  the  conceptual  design.  Tolerance  synthesis  and  analysis  needs  an  exhaustive  functional 
analysis  mechanism  to  make  sure  that  the  identified  functional  requirements  between  the  mating 
components  of  the  assembly  are  met  and  are  suitably  described  in  the  form  of  critical  toleranced 
dimensions  / size  / sizes  / forms  or  in  the  form  of  toleranced  gaps.  The  tolerance  representation 
procedure  uses  the  small  deviation  torsor  scheme  [BAL00]  to  represent  the  variations  associated  with 
each  feature  of  a part  in  the  assembly.  The  tolerance  representation  and  synthesis  scheme  has  been 
elaborated  in  Section  2.  Section  3 discusses  the  formulation  of  cost  functions  while  Section  4 presents 
a discussion  on  the  several  constraint  generation  procedures  for  establishing  the  optimization  process 
as  described  in  section  5.  Section  6 briefly  describes  the  implementation  issues  with  an  example. 


2.  Tolerance  Synthesis  Scheme 

In  this  report  we  propose  the  following  steps  for  the  tolerance  synthesis  scheme: 

1.  For  each  component  of  the  artifact,  develop  the  kinematic  Degrees  Of  Freedom  (DOF) 
requirement  between  the  mating  features  and  decide  the  intrinsic  displacement  and  deviation  of 
the  surface  from  the  nominal  surface  as  torsor  components. 

2.  For  each  pair  of  connected  features,  establish  the  gap  between  the  features  as  a gap  torsor 
taking  into  account  the  indeterminate  and/or  blocked  degrees  of  freedoms. 

3.  Generate  assembleability  constraints  by  considering  all  chains  of  connected  components  as  a 
kinematic  link  (either  closed  loop  or  path  between  two  datum  features)  and  aggregate  the 
torsors  along  the  path  with  reference  to  a global  frame  of  reference.  Each  loop  would  generate  a 
maximum  possible  of  six  equations  binding  the  deviation  parameters  associated  with  all  the 
related  surfaces.  These  sets  of  constraints  would  form  a deviation  for  the  possible  geometric 
variations  of  the  connected  components. 

4.  Generate  functional  requirement  constraints  based  on  the  specified  functions/objectives  to  be 
fulfilled  by  the  artifact.  This  will  restrict  the  range  of  values  of  the  deviations  within  which  a 
component  must  remain  to  fulfill  specific  functional  requirements. 

5.  Generate  manufacturing  / production  cost  models  for  each  component  (based  on  a suitable 
knowledgebase  for  manufacturing/  production  costs)  for  different  machining  operations,  and 
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formulate  a minimization  scheme  by  considering  the  sum  of  costs  for  machining  all  the  parts  as 
the  cost  function  subject  to  the  constraints  defined  in  steps  3 and  4. 

6.  Use  common  mathematical  tools  to  solve  the  above  optimization  scheme  to  find  the  deviation 
parameters  for  each  component.  Convert  the  deviation  parameters  associated  with  each  surface 
into  appropriate  tolerance  zones  through  a tolerance-mapping  scheme. 

Example  for  Illustration 

Let  us  take  an  example  to  elaborate  the  steps  described  above  for  the  tolerance  synthesis  process.  The 
example  has  been  taken  from  the  design  evolution  solution  in  our  earlier  work  [ROYOl]  where  we 
illustrated  the  design  evolution  process  by  taking  an  example  of  designing  a device  for  transmitting 
torque  with  changing  rotational  speed.  We  used  a sample  artifact  library  and  arrived  at  a few  tentative 
design  solutions  that  matched  the  functional  requirements.  We  take  one  of  the  solutions  from  the  above 
design  example  and  assume  the  tentative  sizing  of  the  components. 

The  gear-subassembly  consists  of  three  parts:  a spur  gear,  a shaft  with  a keyhole  and  a key.  The 
geometric  models  of  these  three  parts,  including  their  face-level  mating  relationships  and  behavioral 
models  are  shown  in  Figure  1,  Figure  2. 
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Figure  2:  Mating  surfaces  of  the  gear  sub-assembly 


Part  2 


Figure  3:  Assembly  graph  of  the  gear  sub-assembly 


After  the  torsors  for  each  of  the  mating  surfaces  have  been  defined,  we  consider  the  assembly  graph  of 
the  artifact  (Figure  3)  and  enumerate  closed  loops  of  independent  paths  for  forming  equations  that 
represent  various  possible  geometric  configuration  hulls  for  the  artifact.  Thus,  for  example,  a typical 
closed  loop  from  the  Figure  3 would  be  represented  by:  T[F  1/P  1 ] ->  Gap[Fl/Pl,Fl/P2]  -> 
T[F1/P2,F4/P2]  ->  Gap[F4/P2,F3/P3]  T[F3/P3]  Gap[F3/P3,F4/Pl]  ->  T[F4/P1]  T[F1,P1], 

(see  Appendix-A).  The  aggregation  of  these  torsors  (transforming  into  global  reference  frame,  as 
needed)  should  then  be  equated  to  a null  torsor.  This  will  give  six  equations  connecting  the 
components  of  the  torsors  associated  with  the  loop.  The  gap  between  the  mating  features  will  be  used 
as  constraints  to  satisfy  specific  design  functional  requirements  as  well  as  kinematic  requirements. 
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For  computations  in  each  loop,  the  local  variations  are  to  be  converted  to  a global  coordinate  system  by 
transformation  of  the  corresponding  torsors  in  the  local  frame.  Thus,  in  the  case  of  this  example,  in,  we 
can  take  the  center  of  the  inner  shaft  as  the  reference  origin  and  express  all  the  components  in  terms  of 
vectors  emanating  from  that  point. 

While  traversing  these  loops,  the  gap  between  the  two  mating  features  would  play  a vital  role  in 
deciding  the  type  of  interface  between  the  two  mating  features.  The  gap  torsor  is  required  to  model  the 
kinematic/  functional  specifications  specific  to  the  given  mating  conditions.  However,  we  also  need  to 
assess  the  possibility  of  any  interference  between  the  mating  features  by  applying  the  torsor  on  some 
specific  points  on  the  two  mating  features.  This  would  require  studying  the  torsor  equations  applied  on 
those  strategic  points.  These  strategic  points  could  be  some  inspection  points  to  be  used  in  future 
inspection  of  the  features.  We  have  kept  points  of  kinematic  interest  in  the  artifact  class  definition  so 
that  we  can  create  those  points  as  required  on  specific  surfaces/features.  As  an  example,  in  case  of  a 
rectangular  feature,  we  can  select  the  four  furthest  points  (comer  points)  as  the  desired  points.  For 
circular  features,  we  can  keep  eight  points  uniformly  placed  on  the  circumference  (Figure  4).  It  may  be 
mentioned  that  additional  points  as  desired  may  be  introduced  at  any  stage  to  improve  accuracy; 
however,  the  number  of  equations  and  complexity  of  the  problem  would  increase  with  each  additional 
point. 


Rectangular:  4 points 
(At  the  end  of  each  diagonal) 


Circular:  8 points  at  equal  intervals 


Figure  4:  Points  of  interest  on  features 


After  all  the  loops  have  been  considered  and  corresponding  equations  have  been  formed,  we  will  get  a 
system  of  equations  that  will  define  the  possible  geometric  variations  of  the  artifact  configuration. 
Depending  on  the  type  of  mating  and  functional  requirements  for  each  mating  pairs  of  surfaces,  these 
sets  of  equations  may  be  under-constrained  or  over-constrained.  In  other  words,  there  may  be  too  few 
equations  connecting  too  many  variables  or  vice-versa.  These  will  constitute  the  constraints  of  the 
optimization  process  for  tolerance  synthesis. 

The  proposed  optimization  scheme  could  be  now  summarized  as  below: 
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Minimize:  Cost  of  manufacturing 

Subject  to:  Constraints  associated  with  the  geometry  of  the  product  as  imposed  by  the  kinematic 

requirements  between  mating  features  and  the  constraints  associated  with  functional 
requirements. 

In  the  following  sections,  we  will  detail: 

i)  Proposed  model  for  cost  of  manufacturing 

ii)  Constraints  related  to  Performance/Functional  requirements 

iii)  Constraints  related  to  Assembleability /Kinematic  Mating 

iv)  Constraints  related  to  mapping  of  tolerance  to  feature  deviation  parameters,  and 

v)  Optimization  process,  and  other  implementation  issues 

3.  Generation  of  Constraints 

We  have  considered  two  types  of  constraints  associated  with  the  tolerance  synthesis  process. 

Constraints  arising  to  meet  the  functional  requirements  and  constraints  associated  with  the  kinematic 
mating  requirements  / assembleability.  Since  both  the  constraints  play  a crucial  role  in  the  tolerance 
synthesis  process,  we  have  elaborated  each  process  separately  in  the  following  sub-sections. 

3. 1 Constraints  Related  to  Performance/  Functional  Requirements 

Functional  requirements  are  design  requirements  specified  explicitly  as  well  as  implicitly  by  the  user 
and  designers  to  meet  a certain  desired  performance  from  the  artifact/assembly.  In  terms  of  low-level 
tolerance  and  deviation  parameters  the  functional  requirements  would  translate  into  constraint  i.e., 
functions  establishing  some  relationship  between  parameters.  These  relations  could  be  of  two  types: 
equality  and  inequality. 

The  functional  requirements  can  be  thought  of  in  two  different  ways:  as  a goal  satisfaction  to  minimize 
a penalty  cost  for  deviating  from  the  design  specification  or  to  convert  the  functional  requirements  as  a 
set  of  constraints  on  the  tolerance  zone.  Since  the  functional  requirements  are  developed  based  on 
customer’s  specification  and  designer’s  technical  requirements,  it  is  desirable  to  satisfy  the  product 
specification  / functional  requirements  as  necessary  constraints  rather  than  a goal  satisfaction  process. 
Thus,  we  have  decided  to  treat  the  functional  requirements  as  constraints  on  the  geometric  variations  of 
the  mating  features  of  the  components  and  considered  manufacturing  cost  as  a minimization  goal  to 
arrive  at  an  optimal  tolerance  design.  We  need  to  develop  a generic  process  to  transform  the  functional 
requirements  associated  with  a design  into  a corresponding  set  of  constraints.  The  functional 
requirements  mainly  translate  into  algebraic  inequalities  representing  constraints  of  the  form: 

G(x)  <0,  x eRn.  (3.1-1) 

Functional  requirements  specified  between  two  mating  features  could  be  represented  by  using  the  gap 
torsor  parameters  associated  with  each  mating.  A Functional  requirement  can  be  thought  of  in  two 
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different  ways:  either  as  a goal  satisfaction  to  minimize  a penalty  cost  for  deviating  from  the  design 
specification  or  as  a set  of  constraints  on  the  tolerance  zone. 

We  are  working  to  evolve  a generic  approach  to  transform  the  functional  requirements  associated  with 
a design  into  a corresponding  set  of  constraints.  We  have  already  introduced  some  treatment  of 
functional  requirements  and  constraints  in  our  earlier  work  [ROYOl]  and  we  need  to  extend  further  the 
work  done  to  transform  the  constraint  equations  as  functions  of  the  deviation  parameters. 

3.2  Constraints  related  to  Assembleability  / Kinematic  Mating 

The  constraints  arising  from  the  kinematic  requirements  between  the  mating  surfaces  to  meet  the 
requirement  of  transmission  and/or  blocking  the  desired  flow  of  energy  / force  / fluent  / information 
through  the  mating  area  could  be  treated  as  kinematic  constraints.  The  torsor  scheme  for  the  modeling 
of  the  small  geometric  variations  of  connecting  features  / faces  of  an  artifact  has  been  adopted  for  the 
present  study  and  the  scheme  has  been  elaborated  in  our  earlier  report  [ROYOOa], 

These  kinematic  assembleability  constraints  are  established  by  considering  the  sum  of  the  torsor 
components  along  each  independent  paths/loops  that  exist  between  datum  surfaces.  The  torsor 
summation  is  carried  out  using  the  torsor  transformation  rules  by  transferring  the  effect  to  a global 
point  on  the  datum  surface.  In  our  assumed  model  with  small  displacements/deviations,  the 
assembleability  constraints  become  linear  and  can  be  put  into  the  compact  form:  Ax  = b,  where  A = 
[dij]  an  w x n matrix,  a,y  e R,  x e Rn  b eRm  m is  the  number  of  constraints  and  in  general  m is 
less  than  number  of  deviation  parameters  n.  We  would  need  to  work  further  to  establish  an  automated 
process  to  generates  the  A matrix  from  the  specified  assembly  configuration. 

3.2.1  Derivation  of  the  A matrix 

Since  we  will  the  need  torsor  transformations  in  the  derivation  of  the  A matrix,  the  torsor 
transformations  are  presented  below  (Figure  5). 

Transformational  Torsor 


Figure  5:  Transformation  of  Torsors  - Local  to  Global 
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A small  deviation  torsor  (SDT)  TV  defined  at  a point  /V  in  the  local  coordinate  system  (x,  y,  z)  with  six 
components  ( Qx,  Qy,  Qz,  Sy  , &)  (where  6 denotes  rotational  and  8 denotes  linear  components  of 
movements)  transforms  to  (TN)ML  (where  L denotes  local  coordinate  system)  in  the  same  local 
coordinate  system  shifted  to  a new  point  M is  be  given  by:  (TjQml  = (Ox,  0V,  6Z,  8X  + (Ov*Az  - 6Z  *Az), 
Sy  + (0-_  *Ax  - 6X  *Az),  8:  + (0X  *Ay  - 6V  *Ax))  . This  could  be  written  in  compact  notation  as:  (TjQml  = 
DTn  where  D is  the  6x6  transformation  matrix  given  by 


1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

-A z 

Ay 

1 

0 

0 

A z 

0 

-Ax 

0 

1 

0 

-Ay 

Ax 

0 

0 

0 

1 

(3.2. 1-1) 


where  A =(Ax,  Ay,  Az)r  = (Xn-Xm,  Yn-Ym,  Zn-Zm)t  is  the  vector  from  point  M to  N in  global  coordinate 
system  (X,  Y,  Z). 

After  the  above  transformation,  the  components  of  torsor  (Tn)ml  are  still  in  the  local  coordinate  system 
(x,y,z).  Let  the  local  coordinate  system  (x,y,z)  be  defined  by  the  three  unit  vectors  x = (kxx,  kxy,  kxz), 
y - (kyx,  kyY,kvz)  and  z = (kzx,  kzy,  k2z)  in  the  global  coordinate  system,  in  other  words,  the  local 
coordinate  system  could  be  represented  by  the  [Xj  matrix  given  by: 


M 


(3.2. 1-2) 


The  components  of  (Tn)ml  could  then  be  converted  to  the  global  system  by  the  following 
transformation  (applicable  for  both  the  linear  movement  and  the  small  rotations  in  the  sense  of  SDT) : 


Q x ~ Ox  *kxx+  Qy  *kyx+  Qz  *kzx„ 

(3.2. 1-3) 

Qy  — Qx*  kxy  + Qy*  kyY  + Q*kzy, 

(3.2. 1-4) 

0Z  = Qx*  kxZ  + Qv  * kyZ  + Qz*  kzZ, 

(3.2. 1-5) 

ec  = [0X  eY  ej  = [K][  0, 0, 0JT, 

0g  = [K]  0i. 

(3.2. 1-6) 

where  the  subscripts  g and  l refer  to  global  and  local  coordinate  systems  respectively. 

Sx=  [Sx  + (0y  *Az  - 6z*Ay)]*  kxx+  [Sy+  (0z*Ax  - 6x*Az)]*  kyX+  [&+  (6x*Ay  - 0y*Ax)]*  k^,  (3.2. 1-7) 
8y  — [8X  + (Qy  *Az  - Qz*Ay)]*  kxy+  [8y+  (0z*Ax  - Qx*Az)]*  kyY  + [8Z  + (6x*Ay  - 6y*Ax)]*  kzy  (3.2. 1-8) 
8z=  [8x  + (Qy*Az  - Qz*Ay)]*  kxZ+  [8y+  (6-*Ax  - Qx*Az)]*  kyZ  + [8Z  + (Qx*Ay  - Qy*Ax)]*kzZ  (3.2. 1-9) 

8c,  = [Sx  8y  8zf  = [K]([8x  Sy  SJT  + [Ax  Ay  Azf  x[Qz  Qy  6zf).  (3.2.1-10) 

where  x is  the  vector  cross  product. 
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SG  =[K]  (Sl+  A x 0L)  where  ( 5L=  [S*  5y  <5JT , A =(Ax,  Ay,  Az)T,  0L  = [ 0X  0y  0Z]T ). 

Combining  the  two  transformations,  we  can  write  the  result  as:  (where  subscript  G refers  to  global) 
(Tn)mg  - (Og,Sg) 

(Tn)mg  = ([K][0X  0y  0-f  , [K]  ([8*  5y  S-f  + [Ax  Ay  Azf  x [0X  0y  0-f).  (3.2.1-11) 

The  above  relation  could  be  written  as 

(rj«,  =^2(Tm)m,  (3.2.1-12) 

where  K2  is  the  6 x 6 transformation  matrix  generated  from  (3.2. 1-10). 

We  will  have  to  use  above  results  for  transforming  each  torsor  (defined  in  the  local  coordinate  system) 
in  a loop  to  the  global  reference  point  before  they  can  be  combined.  For  each  loop/cycle  in  the 
assembly  graph  (parts  connected  through  mating  of  features),  we  would  have  (at  most)  six  equations 
by  summing  all  the  deviations  along  a cycle/loop  to  zero  (or  a pre-assigned  quantity).  In  this 
derivation,  we  will  use  the  following  notation.  Without  any  loss  of  generality  and  without  any 
ambiguity  we  will  assume  that  each  loop/path  could  be  enumerated  as  a linked  list  L defined  by: 

£=  {(pj\g^f,pJ\g\-Af’P)}’  (3.2.1-13) 

where  p = part ,/=  feature,  g = gap  between  mating  features.  Symbolically  this  can  be  represented  as 
shown  in  Figure  6. 

Items  within  [ ]...  could  be  repeated  zero  or  more  times  as  required  to  traverse  a loop/path,  (p,f)  is  the 
starting  point  (from_part,  fromfeature),  (f,  p)  the  end  point  (tofeature,  to_part),  [(f,  p,  f),  g]  are 
intermediate  connections  (to  feature,  part,  from  feature),  gap  along  the  path  to  the  next  part. 


Zero  or  more  reDetitions 


Part 

Feature 

✓ 

* 

/ 

Gap  { 

Feature  Part  Feature 

\ 

\ 

Gap  ; 

Feature  Part 

s 

✓ 

Figure  6:  Loops/Paths  in  an  assembly 


In  case  of  a closed  loop,  the  last  part  would  be  the  same  as  the  first  part.  Both  the  starting  and  the  final 
parts  for  open  loops  should  be  terminals  (datums).  In  the  case  of  closed  loops,  the  summation  of  the 
torsors  would  be  null  (0),  whereas  for  an  open  path,  the  summation  would  be  equal  to  a torsor 
(possibly  non-null)  that  has  to  be  established  from  known  functional  specifications  between  the  two 
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terminal  parts.  We  will  designate  this  quantity  as  b.  We  will  use  T to  indicate  torsor  with  index  p,fg 
for  parts,  features  and  gaps  respectively. 

Pj  = Part  number  j,  j = 1 to  n,  n-  number  of  parts  in  the  loop, 

Fk  = Feature  k of  part  Pj , k = 1 to  m,  m = number  of  features  of  part  Pjt 

Tp  = 6-component  displacement  torsor  of  the  part  P,  defined  at  a suitable  point  O inside  the  part  in  a 
local  coordinate  system, 

Tf  = 6-component  deviation  torsor  of  the  feature  k (of  some  part  Pj) , defined  at  a suitable  point  N on 
the  feature  in  a local  coordinate  system  with  x-axis  along  the  outward  (away  from  the  material) 
normal  and  y and  z are  the  other  two  axes  mutually  perpendicular  to  the  x-axis, 

Tg  = 6-component  gap  torsor  defined  at  a mating  between  two  features  of  two  parts. 

With  the  above  notation,  we  can  write  the  desired  equations  for  each  loop,  with  implied  local  to  global 
transformations  applied  to  each  torsor  as  viewed  from  any  convenient  point  (typically,  at  the  center  of 
the  first  part  where  the  loop  begins),  as: 


(Tp+  Tf)  + Tg  /+  (Tf  + TP  + Tf  + Tg)]...  + (Tf+  Tp)  = b. 


(3.2.1-14) 


Summing  over  the  entire  loop,  we  can  write, 

J^Kj*(Sj+Axjx0j)  = bs, 

j e EntireLoop 

j e EntireLoop 


(3.2.1-15) 


(3.2.1-16) 


where, 

Kj  is  the  local  to  global  transformation  matrix  (eqn.  (3.2. 1-2) ) for  feature  j, 

Sj  = [Sx  5y  Sz]Tj  and  0j  = [9z  9V  0z]Tj  are  the  two  vectors  formed  from  the  linear  and 
rotational  components  respectively  of  the  deviation  torsor  of  the  feature  j. 

The  scheme  for  the  representations  and  transformations  of  the  torsors  have  been  elaborated  in 
Appendix  - A with  an  example  of  an  artifact  with  three  blocks. 

Since  the  proposed  scheme  uses  six  small  displacement  parameters  to  represent  the  possible  geometric 
deviations  from  the  design  / nominal  surface,  we  need  a mechanism  to  express  standard  tolerances  in 
terms  of  these  torsor  parameters.  We  consider  some  engineering  features/surfaces  and  develop  the 
tolerance  to  deviation  mappings  in  the  following  section. 
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3.3  Constraints  related  to  mapping  of  tolerance  to  deviation  parameters 

The  category  of  relations  arising  from  the  mapping  of  the  deviation  parameters  of  the  features  to 
corresponding  tolerance  values  is  of  the  generic  non-linear  inequality  (“less  than  or  equal  to”  type)  and 
takes  the  form:  Fj>p  (x)  <FTok  (t),  where  FDP  is  some  function  of  the  deviation  parameters  x eR",  Ft0i 
is  some  function  of  the  tolerance  parameters  t , and  ke  N*  is  an  index.  Exact  forms  of  Fdp(x)  and 
Froi(t),  depend  on  the  nature  of  the  feature/surface.  For  a planar  surface,  the  above  relationship 
becomes  a set  of  linear  constraints.  The  interesting  feature  in  the  above  relationship  is  the  separable 
nature  of  the  equation  in  terms  of  the  variables  t and  x.  We  would  like  to  explore  the  possibility  of 
taking  advantage  of  the  separable  nature  of  the  mapping  to  define  the  deviation  hull  in  a semi- 
analytical  way,  if  possible.  The  linear  nature  of  the  equations,  however,  does  not  exist  when  we 
consider  other  non-planar  features.  So  far,  we  have  been  successful  in  establishing  the  required 
relationships  for  the  features:  rectangular  planar,  circular  planar,  cylindrical,  and  spherical.  The 
following  sections  give  the  detailed  results. 

3.3. 1 Representation  of  Geometric  Deviations  of  Features  of  a Part 

The  geometric  variations  of  the  features  are  defined  in  terms  of  deviations  of  their  independent  degrees 
of  freedom.  In  this  approach,  at  the  most  six  independent  parameters  would  be  used  to  represent  a six- 
component  vector  corresponding  to  the  possible  six  degrees  of  freedom.  In  actual  case,  depending  on 
the  type  of  surface  (planar,  cylindrical,  etc)  some  of  the  degrees  of  freedom  are  ‘free’  (variations  along 
which  keeps  the  surface  invariant)  and  thus,  in  general,  less  than  six  deviation  parameters  are  required. 
For  example,  for  a planar  surface,  rotation  along  the  outward  normal  to  the  surface  and  movements 
along  the  two  axes  on  the  surface  keep  the  feature  invariant  and  hence  a planar  surface  requires  three 
deviation  parameters  (axial  movement,  two  rotations  along  the  two  transverse  axes). 

In  order  to  make  the  process  of  converting  the  functional  requirements  as  well  as  kinematic 
requirements  into  a set  of  constraints  in  a generic  manner,  we  need  a procedure  to  represent  the 
variations  associated  with  each  part  in  terms  of  a set  of  generic  parameters  of  the  mating 
features/surfaces.  Assuming  small  displacements/deviations  of  these  features  due  to  manufacturing 
inaccuracies/other  defects,  the  actual  surface  of  a feature  deviates  from  the  theoretical  (nominal) 
surface  slightly.  Such  variations  could  be  systematically  treated  by  considering  variations  of  six  small 
parameters  defined  at  a point  on  the  surface.  The  six  parameters  are:  three  linear  translations  and  three 
rotations  along  the  three  axes  of  the  local  orthogonal  coordinate  system  corresponding  to  the  six 
degrees  of  freedom.  These  six  small  displacement  parameters  have  been  called  screw  parameters  or 
torsors. 

Since  the  proposed  method  uses  six  small  displacement  variation  parameters  to  represent  the  possible 
geometric  deviations  from  the  design  / nominal  surface,  we  need  a mechanism  to  express  geometric 
variations  in  terms  of  these  torsor  parameters.  It  has  been  mentioned  above  that  not  all  the  six 
parameters  for  each  surface  may  be  uniquely  defined  and  we  may  not  need  to  know  all  the  six 
parameters  in  specific  cases.  We  plan  to  consider  each  type  of  engineering  surfaces  and  formulate 
methods  to  represent  the  tolerance  values  as  functions  of  the  parameters.  This  task  require  a mapping 
of  the  deviation  parameters  to  geometric  tolerances.  So  far,  we  have  completed  mappings  for  Planar, 
Cylindrical,  and  Spherical  features.  Conical,  Toroidal,  and  Helical  cases  will  be  taken  as  future  work. 
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In  order  to  give  a meaningful  interpretation  to  the  assemblability  of  an  artifact  in  terms  of  tolerance 
specifications  as  per  industry  standard  codes  and  practices  like  American  Society  of  Mechanical 
Engineers  (ASME)  Y14.5.1M-1994  [ASMEY14],  International  Organization  for  Standardization  ISO 
[ISO],  we  need  to  consider  relationship  between  the  tolerance  zones  (or  virtual  condition  boundaries) 
and  the  deviation  parameters  of  a feature.  By  deviation  parameters,  we  mean  the  representation  of  the 
possible  deviation  of  a feature  by  using  (at  most)  six  independent  parameters  corresponding  to  the  six 
degrees  of  freedom  of  a feature.  In  general,  the  deviation  parameters  to  tolerance  mappings  are  not 
one-to-one  and  can  only  be  represented  as  inequalities  involving  functions  of  the  tolerance  parameters 
and  the  deviation  parameters.  Mappings  developed  from  these  are,  in  general,  a set  of  non-linear 
inequalities.  However,  for  some  specific  cases,  like  the  planar  surface,  the  inequalities  are  linear,  and 
represent  a diamond  shaped  zone  by  the  intersection  of  inequalities,  each  of  which  are  half-spaces  in 
the  D x T ( D : Deviation,  T:  Tolerance)  space  [BAL98]. 

There  are  four  types  of  tolerances  used  in  the  industry:  size  tolerance,  form  tolerance,  orientation 
tolerance,  and  location/position  tolerance.  These  tolerances  are  sometimes  specified  with  material 
conditions  as  modifiers.  There  are  three  such  material  conditions:  MMC  (Maximum  Material 
Condition),  LMC  (Least  Material  Condition)  and  RFS  (Regardless  of  Feature’s  Size).  Since  we  assume 
that  the  shape  of  a toleranced  feature  remains  similar  to  the  nominal  feature  (for  example,  a plane 
remains  a plane,  a cylinder  remains  a cylinder,  etc.),  the  variation  of  position  and  orientation  due  to 
form  tolerance  is  very  small  [Gil91]  and  in  this  study,  we  consider  only  the  size  tolerance,  orientation 
tolerance,  and  positional  tolerance. 

A geometric  tolerance  applied  to  a feature  of  size  and  modified  to  MMC  will  establish  a virtual 
condition  boundary  (VCB)  4 outside  of  the  material  space  adjacent  to  the  feature.  The  feature  shall  not 
cross  this  VCB.  Likewise,  a geometric  tolerance  applied  to  a feature  of  size  and  modified  to  LMC 
establishes  a VCB  inside  the  material  and  the  feature  shall  not  cross  the  VCB.  This  translates  to  the 
following  relationships  [DRA99]: 


Modified  to  MMC: 

For  an  internal  feature  of  size: 

MMC  virtual  condition  = MMC  size  limit  - geometric  tolerance;  (3.3. 1-1) 

For  an  external  feature  of  size: 

MMC  virtual  condition  = MMC  size  limit  + geometric  tolerance;  (3.3. 1-2) 

Modified  to  LMC: 

For  an  internal  feature  of  size: 

LMC  virtual  condition  = LMC  size  limit  + geometric  tolerance;  (3.3. 1-3) 

For  an  external  feature  of  size: 

LMC  virtual  condition  = LMC  size  limit  - geometric  tolerance.  (3.3. 1-4) 

If  a geometric  tolerance  applied  to  a feature  of  size  is  modified  to  neither  MMC  nor  LMC,  by 
definition  of  standards  such  as  ASME  Y14.5,  it  is  modified  to  RFS.  In  this  case,  instead  of  VCB,  the 
tolerance  specification  will  generate  a tolerance  zone  into  which  the  derived  element  will  not  interfere. 
A VCB  or  tolerance  zone  represents  the  basic  intention  of  the  designer  that  each  point  on  the 
toleranced  feature/derived  element  should  remain  within  (or  outside)  this  boundary/zone.  Hence,  the 
VCB  or  tolerance  zone  would  be  ideally  suitable  for  representing  the  relational  limits  on  the  deviation 


4 Virtual  condition  boundary:  A constant  boundary  generated  by  the  collective  effects  of  a size  feature’s  specified  MMC 
or  LMC  material  condition  and  the  geometric  tolerance  for  that  material  condition. 
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parameters  associated  with  the  deviation  of  the  feature  from  its  nominal  position.  Since  the  deviation 
parameters  for  a feature  could  be  used  to  define  the  deviation  of  all  points  on  a feature  (by  employing 
one  or  more  independent  parameters),  the  tolerance  specifications  could  be  thought  of  as  a set  of  limits 
for  the  deviation  parameters  and  vice-versa.  These  relationships  could  be  used  for  tolerance  synthesis 
as  well  as  for  tolerance  analysis  including  checking  the  assemblability  in  the  worst-case  scenario. 

In  the  following  sections  of  this  paper,  we  consider  the  mapping  for  size  tolerance,  positional 
tolerance,  and  orientation  tolerance  for  four  types  of  features:  planar,  cylindrical,  spherical,  and 
conical.  For  the  planar  case,  there  is  no  material  condition  as  it  is  a non-size  feature.  For  the  cylindrical 
feature  we  have  two  cases:  MMC  and  RFS  and  for  the  spherical  feature  we  have  MMC  and  for  conical 
feature  we  have  the  size  tolerance  only. 

We  use  the  following  basic  steps  to  convert  a tolerance  specification  into  a set  of  inequalities  in  the 
deviation  parameters: 

a.  Generate  the  intrinsic  torsor  and  the  deviation  torsor  for  the  feature  by  eliminating  the  deviation 
parameters  that  are  invariant  for  the  feature  to  reduce  the  degrees  of  freedom.  This  reduces  the 
number  of  deviation  parameters  needed  to  represent  the  variation  of  the  feature.  For  example, 
for  a cylindrical  feature,  axial  movement  and  axial  rotation  are  two  invariants.  Hence,  we  do 
not  need  these  two  parameters. 

b.  Generate  a VCB  (or  a tolerance  zone)  based  on  the  tolerance  specification.  In  this  step,  we 
compute  the  size  of  the  VCB  or  the  tolerance  zone  for  restricting  the  variation  of  the  feature  or 
derived  element  respectively. 

c.  In  case  of  VCB,  take  an  arbitrary  point  on  its  nominal  surface  of  the  feature  in  a parametric 
form  and  transform  it  to  a new  position  by  applying  the  effect  of  the  deviation  torsors.  For 
example,  for  a cylindrical  feature  a point  P on  its  nominal  surface  is  represented  by  two 
parameters  (9,  z);  P = (rcosO , rsinO,  z),  where  6e(0,  2ir)  and  z e (0,  L),  L=Length  of  the 
cylinder.  In  the  case  of  tolerance  zones,  transform  the  derived  element,  such  as  center  plane, 
center  axis  etc.,  to  a new  position  by  applying  the  effect  of  the  deviation  torsor. 

d.  Eliminate  the  free  parameters,  (such  as  the  (6,  z)  mentioned  in  step  c above),  by  applying  the 
condition  that  the  extreme  points  of  the  transformed  position  should  remain  within  the  VCB  or 
that  the  derived  element  should  remain  within  the  tolerance  zone.  This  step  generates  a set  of 
inequalities  connecting  the  deviation  parameters  with  the  tolerance  specification. 

The  procedure  for  applying  the  above  four  steps  will  be  clear  in  subsequent  sections  where  we  generate 
the  mappings  for  the  basic  features  under  various  material  conditions. 

The  symbols  used  in  the  following  sections,  unless  otherwise  specified,  are  defined  as  below: 

TU  = Upper  limit,  size  tolerance, 

TL  = Lower  limit,  size  tolerance, 

Tp  = Positional  tolerance, 

Tv-  Perpendicularity  tolerance, 

(A9x,  A9y,  A0Z,  Ax,  Ay,  Az)  = six  components  of  the  deviation  torsor. 
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3.3.2  Planar  Features 

For  each  surface/feature,  a local  coordinate  system  (LCS)  is  defined  and  the  deviation  parameters  are 
defined  in  that  LCS.  For  a planar  surface,  the  LCS  is:  the  z-axis  is  the  outward  normal  (emanating  from 
the  material  side  of  the  feature)  and  (x,  y)  are  local  orthogonal  coordinates  on  the  plane,  so  that  the 
equation  for  the  nominal  surface  (plane)  is  given  by  z - 0.  For  this  plane  the  deviation  parameters  of 
the  small  deviation  torsor  (SDT)  are:  (A Qx,  A6y,  0,  0,  0,  Az)T. 

Size  Tolerance  1:  Size  tolerance  with  respect  to  a datum. 

Case  1:  Rectangular  (2a  x 2b)  planar  surface.  (Figure  7) 


Figure  7:  Planar  feature  (rectangular)  with  size  tolerance 


For  the  above  feature,  deviation  parameters  are  (A6Xi  A0y  0,  0,  0,  Az),  and  the  tolerance  parameters  are 
(TU,  TL)  and  the  SDT  is:  Dd  = {AO  Ad}7  = {A6X  A6y  0 0 0 Az}7  (3.3.2-1) 

Since  this  is  a rectangular  planar  feature,  the  four  comers  with  nominal  coordinates  (±a,  ±b,  0)  will 
have  the  maximum  deviation  due  to  the  effect  if  Dd.  Transforming  Dd  on  to  these  four  points  we  get, 
D ’ = Dd+  ((±a,  ±b,  0)  * (A6X  A6y  0)  which  gives  deviations  for  the  four  points  as: 

D,  = {A0X  AfyO  0 0 Az  + a*A0y- b*A0x }T, 

D2  - { A6X  A0V  0 0 0 Az  + a*A6y  + b*A6x }7, 

D3  = {A0x  A9y  0 0 0 Az-  a*A0y+  b*A6x }T, 

D4  = {A0X  A0yOO  0 Az-  a*A6y  - b*AGx  }r. 

Comparing  the  z-component  with  the  specified  tolerance  value  (TU,  TU),  we  get  the  following  four 
inequalities  for  the  desired  mapping: 

-TL  < Az  + a*A6y-  b*A0x  < TU,  (3. 3.2-2) 
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-TL  < A + a*A0v  + b*A9x  < TU, 

-TL  < A-  a*A0v  + b*A0x  < TU,  and 
-TL  < A - a*A0y  - b*A0x  < TU. 


(3. 3.2-3) 
(3. 3. 2-4) 
(3. 3. 2-5) 


Case  2:  Circular  (radius  r ) planar  surface.  (Figure  8) 


Figure  8:  Planar  feature  (circle)  with  size  tolerance 


Deviation  parameters  = (0,  0,  A,  A6X  AGy  0),  tolerance  parameters  = (TU,  TL). 

For  arbitrary  point  P on  the  circumference,  at  an  angle  / 3 with  the  y-axis,  we  have, 

Sz  = A -RA6vSin(/3)+RA0xCos(J3), 

where  Sz  is  the  pure  position  variation  of  a point  in  the  Z direction,  where  there  is  a size  tolerance 
control.  The  tolerance  zone  in  this  case  is  made  up  of  two  planes, 

z = +TU  and  z = -TL. 

We  need  to  ensure  that  4 is  within  the  above  boundary  V(3  e (0,  2n) 

This  can  be  done  by  treating  the  8Z  as  a function  of  (3  and  then  finding  the  extreme  points  by  putting 
3(8- )/  dp  = 0.  This  gives, 

8 max/mm  A ± R V(  A0X"  +A0y~). 


This  leads  to  the  following  mapping  (constraints): 

(A0X  + A9y2)  < (TU+  TL)2/(4R2). 
We  can  also  write  in  the  parametric  form  as 
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TL  < Az  + RA6xSin(p)  + RA0yCos(/3)  <TU , V(3  e ( 0 , 2n). 

The  last  inequality,  using  the  value  of  Sz  max/mm  could  be  written,  independent  of  /?,  as: 

TL<  (Az  -R  V(  A6X2  +A0y2)  ) and  (Az  + R V(  A 6X  +A9y2)  ) < TU . 

3.3.3  Cylindrical  Feature  - Modified  to  MMC 

Let  us  assume  a cylindrical  feature  with  a given  tolerance  specification  as  shown  in  Error!  Reference 
source  not  found..  A local  coordinate  system  (LCS),  as  shown  in  Error!  Reference  source  not 
found.,  is  defined  for  cylindrical  surface  and  the  deviation  parameters  are  defined  in  that  LCS.  For  a 
cylindrical  surface,  the  LCS  is:  z-axis  along  the  axis  of  the  cylinder  and  (x,  y)  are  local  orthogonal  co- 
ordinates on  the  middle  of  the  axis  (Error!  Reference  source  not  found.,  Figure  10).  For  this  surface, 
the  deviation  parameters  of  the  SDT  are:  (Ax,  Ay,  0,  A0X,  A6y,  0)  and  dr.  Based  on  this  notation, 
constraints  could  be  derived  connecting  these  parameters  with  the  specified  tolerance  values  as 
detailed  below. 

Deviation  parameters  = (Ax,  Ay,  0,  AGXi  A By,  0),  tolerance  parameters  = (TU,  TL,  Tp),  SDT  is: 

Dd  = {AG  Ad}  = { AGX  AGy  0 Ax  Ay  0}T , 

D,  = {Ad  Ad}  ={0  0 0 dr*cos6  dr*sinO  0}T , 

where  Dd  is  displacement  torsor 5 and  D,  is  intrinsic  torsor6. 

Let  us  assume  that  there  is  a point  P (rcosO,  rsinO,z)  (Error!  Reference  source  not  found.)  on  the 
nominal  cylindrical  surface.  After  applying  the  transformation  due  to  the  effect  of  the  two  small 
displacement  torsors,  the  new  position  of  transformed  point  P ' can  be  calculated  as  follows. 


5 Displacement  torsor:  torsor  for  the  deviation  of  the  feature  from  the  nominal  position. 

6 Intrinsic  torsor:  torsor  to  represent  the  intrinsic  variation  of  size  of  the  feature. 
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1 

0 

A*v 

Ax 

(r  + dr) cos  0 

(r  + dr)  cos  0 + A 6yz  + Ax 

0 

1 

-A*x 

Ay 

(r  + dr)  sin  0 

(r  + dr)  sin  6 - A 6xz  + Ay 

-a  ey 

A 9X 

1 

0 

z 

z — A0y  (r  + dr)  cos  0 + A0X  (r  + dr)  sin  0 

0 

0 

0 

1 

1 

1 

The  4 by  4 matrix  on  left  hand  side  of  the  equation  is  the  transformation  matrix  for  point  P following 
the  assumed  transformation  sequence  Ax ->Ay ->A6x->A6y  by  displacement  torsor.  The  4 by  1 matrix 
(vector)  on  the  left  hand  side  of  the  equation  represent  the  changed  position  of  point  P by  intrinsic 
torsor. 


So,  the  new  position  P ' (Error!  Reference  source  not  found.)  is  given  by 
((r+dr)cosj9+A6y  z+Ax,(r+dr)smO-A6x  z+Ay,z-A6y  (r+dr)cos0+A0x(r+dr)smO). 

The  VCB  of  the  cylindrical  surface  is  constructed  by  using  the  rules  (3. 3. 1-1,  3.3.1  -2,3 .3. l-3,3.3. 1 -4) 
and  has  the  following  properties: 


o A perfect  shape  as  that  of  the  nominal  cylindrical  feature, 

o A size  (diameter)  of  maximum  size  plus  geometric  tolerance,  because  its  external  feature  and 
the  positional  tolerance  is  modified  to  MMC, 
o A perfect  orientation  as  that  of  the  nominal  cylindrical  feature  (vertical  to  C), 
o A perfect  location  as  that  of  the  nominal  cylindrical  feature. 


Hence,  the  constraint  equation  will  be: 

((r  + dr) cos#  + A 6yz  + Ax)2  + ((r  + dr) sin  0 - AO xz  + Ay)2  < (2r  + TU  +Tp)2  /4 , (3. 3. 3-1) 

where  6 e (0,  2n).  The  left  hand  side  is  a function  of  two  independent  parameters  6 and  z e (0,  L).  The 
inequality  should  be  valid  V6  e (0,  2n)  and  Vz  e (0,  L).  We  thus  need  to  eliminate  Onnd  z from  the 
LHS  of  (3. 3.3-1)  by  finding  the  maximum  of  the  above  equation.  LHS  of  equation  (3. 3. 3-1)  could  be 
written  as: 

f(0,z)  = {a cos 6 + bz  + c)2  + (asinO  + dz  + e)2 , (3.3.3  - 2) 

where  a = r + dr,  b = A0y,  c-  Ax,  d = -A0X,  e-Ay. 


By  equating  the  partial  derivatives  of  f to  zero, 

— = 0,  and  — = 0,  we  get, 

80  dz 

r ±1  2 2 ^ 

(< ad 2 + ab~)  - be  + cd 


= max< 


yjd2  +b: 


d2  + b2 


,0 


, f - min 

5 J min 


±1 


Jd2+b: 


(ad 2 + ab2 ) - be  + cd 


d2  +b- 


,0 


(3.3.3 -3) 
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Figure  10:  Cylindrical  Feature  with  Tolerance  Modified  to  MMC 


The  final  set  of  constraints  is: 


f < 

J max 


2r  + TU  + Tp 


-TL<  dr  < TU . 

When  the  toleranced  feature  is  internal,  as  for  example,  for  a hole,  the  constraints  are: 

r2r  + TU-Tp' 


f > 

J min 


-TL<  dr  < TU . 


J 


If  the  external  feature  is  modified  to  LMC,  the  constraint  will  be: 


f > 

J mm 


2r-TL-Tp 


-TL<  dr  < TU . 

If  the  internal  feature  is  modified  to  LMC,  the  constraint  will  be: 
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(2  r-TL  + Tp^ 

J max  — ’ 

V 2 ) 

-TL<  dr  < TU . 


3.3.4  Cylindrical  Feature  Modified  to  RFS 


Figure  11:  Cylindrical  Feature  with  Tolerance  Modified  to  RFS 

i 

i 

XI  i 

t L 


Figure  12:  Cylindrical  tolerance  zone 
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Figure  12  shows  a cylindrical  tolerance  zone  (the  cylindrical  feature’s  center  axis)  that  could  be  due  to 
positional  tolerance,  concentricity  tolerance,  and/or  run-out  tolerance  for  the  cylindrical  surface.  A 
LCS  (local  coordinate  system)  is  defined  as  shown  in  Error!  Reference  source  not  found..  The 
displacement  torsor  and  the  intrinsic  torsor  are: 

Dd  = {AO  Ad}  = {AQX  AO,  0 Ax  Ay  0}T, 

D,  = {Ad  Ad}  = {0  0 0 dr*cosd  dr*sind  0}T. 

The  basic  difference  between  this  case  and  the  one  in  Figure  10  is  that  Dd  is  independent  from  D,.  It 
means  that  the  size  of  the  tolerance  zone,  r,  will  not  change  due  to  the  variation  of  dr,  which  is  the 
departure  from  nominal  radius.  Because  it  is  the  derived  element,  axis  of  the  feature  of  size,  that  is 
toleranced  but  not  the  surface  of  the  feature  of  size,  different  methods  are  used  to  generate  the 
constraints  on  mapping  deviation  parameters  to  the  tolerance  zone.  Since  Dd  and  D,  are  independent, 
we  need  to  treat  them  separately.  For  Dd  we  have  4 deviation  parameters:  Adx,  Ady  Ax,  and  Ay.  For 
ease  of  computation,  let  us  assume  that  the  deviations  of  the  axis  takes  place  in  the  sequence  of 
Ax  ->Ay  ->Adx->Adv.  We  think  this  assumption  will  not  affect  the  generality  because  it  does  not  matter 
what  sequence  we  choose.  (In  fact,  the  matrix  multiplication  described  earlier  is  based  on  the  same 
assumption.) 

The  constraints  on  Ax  and  Ay  are: 


- r < Ax  < r 

- Jr2  - Ax2  < Ay  < * Jr 2 - Ax2 


(3. 3.4-1) 


where  r is  the  radius  of  tolerance  zone. 


Figure  13  shows  the  movement  of  the  center  axis  by  Ax  and  Ay  followed  by  a rotation  of  the  center 
axis  by  Adx  until  the  end  of  it  touches  tolerance  zone  boundary. 

We  determine  the  constraint  on  Adx,  given  those  of  Ax,  Ay,  as  follows,  following  Figure  13: 


Figure  13:  Derivation  of  constraint  on  Dqx 


■Jr2  - Ax2  - M <A0  ^r2~A*2  ~N 


L/2 


Lt  2 


(3. 3.4-2) 


Figure  14  shows  the  movement  of  the  center  axis  by  Ax  and  Ay  followed  by  a rotation  of  the  center 
axis  by  Adx,  then  rotate  the  center  axis  by  Adv  until  the  end  of  it  touches  tolerance  zone  boundary. 
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From  Figure  14,  we  get. 


^-(A6>-Z/2+|Aj])2-|Ax|  Jr2  -(At?  -L/2+|Avj)2  — |Axj 


(I/2)cosA^ 


<A 0 <- 

v - 


(L/2)cosAQ 


Because  zl#  is  very  small,  cos(AOx)  &1,  so  we  have, 


^2-(A^-JL/2+|4j)2-|Axj  ^2-(A^L/2+|Av])2  -|Ax| 

^ Ar/  ^ 


(1/2) 


a/2) 


(3. 3.4-3) 


Figure  14:  Derivation  of  constraint  on  Dqy 

Inequalities  (3. 3.4-1, 3. 3. 4-2,3. 3.4-3)  are  the  mapping  results  for  the  positional  tolerance  modified  by 
RFS  specified  on  a cylindrical  feature  of  size. 


3.3.5  Spherical  Surface 

A local  coordinate  system  (LCS)  is  defined  for  a spherical  surface  and  the  deviation  parameters  are 
defined  in  that  LCS.  For  a spherical  surface,  the  LCS  is:  the  z-axis  is  along  a radius  of  the  sphere,  and 
(x,  y)  are  local  orthogonal  co-ordinates  at  the  center  of  the  sphere  (Figure  15,  Figure  16).  For  this 
surface  the  deviation  parameters  of  the  SDT  are:  (Ax,  Ay,  Az,  0,  0,  0 ) and  dr.  Based  on  this  notation, 
constraints  could  be  derived  connecting  these  parameters  with  the  specified  tolerance  values  as 
detailed  below. 


Page  22 


Tolerance  Synthesis  Scheme 


Page  23  of  60 


Figure  15:  Spherical  feature  with  LCS 


S <J>  2 r ±TL 


S0TP® 


Figure  16:  Spherical  feature  with  tolerance  specification 


The  VCB  of  the  spherical  surface  is  constructed  with  the  following  properties: 
o A perfect  shape  as  that  of  the  nominal  spherical  feature, 

o A size  (diameter)  of  2r+TU+Tp,  because  it  is  an  external  feature  and  the  positional  tolerance  is 
modified  to  MMC, 

o A perfect  location  as  that  of  the  nominal  cylindrical  feature. 

By  applying  the  same  procedures  as  was  followed  for  the  cylindrical  case,  we  have: 

Dj  = {AO  Ad}  = {0  0 0 Ax  Ay  Az}7 

D,=  {AO  Ad}  = {0  0 0 dr*sin(j)*cos  y/  dr*sin<f)*sin y/  dr*cos(/>}T,  and 
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'1 

0 

0 

Ax' 

(r  + dr)  sin  $ cos  y/ 

Ax  + (r  + dr)  sin  (j)  cos  y/ 

0 

1 

0 

Ay 

(r  + dr)  sin  (j)  sin  y/ 

Ay  + (r  + dr)  sin  (j)  sin  y/ 

0 

0 

1 

Az 

(r  + dr)  cosy/ 

Az  + (r  + dr)  cos  y/ 

0 

0 

0 

1 

1 

1 

The  4 by  4 matrix  on  the  left  hand  side  of  the  above  equation  is  the  transformation  matrix  for  any  point 
on  the  spherical  surface  due  to  displacement  torsor  assuming  the  transformation  sequence 
Ax  ->Ay  ->Az.  The  4 by  1 matrix  (vector)  on  the  left  hand  side  of  the  equation  represents  the  changed 
position  of  the  point  due  to  the  intrinsic  torsor. 

The  constraint  equation  is  given  by: 

(Ar  + (r  + a,r)sin^cos^/)2  + (Ay  + (r  + dr)sin^sin^)2  + (Az  + (r  + dr)cosy/)2  < (2 r + TU  + TP )2  /4 
and  dr  > -TL  . 

We  need  the  above  inequality  to  be  valid  V0  e (0,  2n)  and  VP  e (0,  2n).  So  we  eliminate  parameters 
0and  Pby  maximizing  the  following  function. 

y/)  - a sin2  (j)  + a cos2  y/  + 2 Ax  sin  ^ cos  y/  + 2 A y sin  (j)  sin  y/  + 2 Az  cos  y/ , 

where  a = r + dr 

df  / d<fi  = 0 2 cos  ^(<2  sin  ^ + Axcos^/  + Aysin^)  = 0, 

dfldy/- 0 =>  - 2<?  cos  ^ sin  ^ - 2 Ax  sin  ^ sin  ^/ + 2 Ay  sin  ^ cos  - 2 Az  sin  = 0. 

Here,  we  have  two  sets  of  solutions  (a  symbolic  software  package  was  used  to  get  the  solutions  in 
symbolic  form): 


1) 


* = (2n  + l)|. 


¥ = ~tg 


a(3  - Ay 
Ax  + Az  ’ 


where  (52  = (z4cr2  - 2z2  a Ay  + (Az2  + Ay2  + 2AxAz  -a2  + A x2)z2  + 2azAy  - Ay2) 


' Ax(-y2 a2  + y2Ax2  + y2Ay2  -yAza-Ay2)' 

( AyAx(2y2  - 1)  \ 

bill 

v (-ya2  + yAx2  - yAy2  - Aza)a  y 

~lg 

v-  y2a2  + y2 Ax2  -y2 Ay2  - yAza ^ 

2)  where  y2  = ((<24  — 2a1  Ax2  + 2a2 Ay2  + Ax4  + Ay4  + 2Ay2Ax2)z 4 + (2A zo3  - 2aAzAx2  + 

2aAzAy2)z 3 + (- a 4 - 2Ay2Ax2  - 2a2 Ay2  - Ay4  - Ax4  + 2a2  Ax2  + A z2a2)z2  + ( 2aAzAx 2 - 2A za2  - 

2aAzAy2)z  - Az2 a2  +Ay2Ax2). 

The  final  set  of  constraints  is: 

f2r  + 7T/  + 7>V 

J max  — ’ 

V 4 J 
-TL<  dr  < TU . 
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When  the  toleranced  feature  is  internal,  for  example  a hole,  the  final  constraint  will  be: 

f 2r  + TU-TpV 


f > 

J min 


v 


-TL<  dr  < TU . 

If  the  external  feature  is  modified  to  LMC,  the  constraint  will  be: 


f > 

J min 


2 r-TL-  Tp 


- TL<  dr  < TU. 

If  the  internal  feature  is  modified  to  LMC,  the  constraint  will  be: 

, (2r-TL  + Tp~ 

f < 

J max  rx 

V ^ 

-TL<  dr  < TU . 


3.3.6  Conical  Feature 

A local  coordinate  system  (LCS)  is  defined  for  a conical  surface  and  the  deviation  parameters  are 
defined  in  that  LCS.  For  a conical  surface,  the  LCS  is:  the  z-axis  is  along  the  cone  axis  and  (x,  y)  are 
local  orthogonal  co-ordinates  on  the  bottom  plane  of  the  cone  (Figure  17).  For  this  surface  the 
deviation  parameters  of  the  SDT  are:  (Ax,  Ay,  Az,  A6X,  AdVi  0)  and  dr.  Based  on  this  notation, 
constraints  could  be  derived  connecting  these  parameters  with  the  specified  tolerance  values  as 
detailed  below. 


Figure  17:  Conical  feature  with  LCS 
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A 0X  Ax 

0 

dr  cos  0 

A6y  Ay 

> D,  = ■ 

0 

dr  sin  6 

0 Az 

0 

0 

(3. 3.6-1) 


A point  P on  the  nominal  surface 


R -r  ^ 

r cos  0,  r sin  6, — ! L 

R,-R2  , 


, will  transform  to: 


1 


0 AG„  Ax' 


0 1 - A Gx  Ay 

- A 6 AG,  1 A z 


0 


0 


0 


1 


(r  + dr)  cos  0 
(r  + dr)  sin  6 

R'~rL 


R\  R2 


(r  + dr)cosG  + AG  — — — L + Ax 
yR,-R2 

(r  + dr)sinG  - Ad 


R\  -R2 


L + Ay 


R -r 

— ! L - A0v(r  + dr)cos6  + A 0x(r  + dr)sind  + Az 

R\~Ri 


For  the  tolerance  specification  of  the  following  feature,  the  constraints  will  be: 


(3. 3. 6-2) 


r-TL<  V*  2 + y'2  + Rl--2  z'<r  + TU , 

L 

x'  - (r  + dr)cos0  + A6  — — —L  + Ax, 
yRx-R2 
R — r 

y - (r  + dr)sind  - A6X  — '■ L + Ay, 

where  Ri~R2 

R -r 

z = — ! L - A 0 (r  + dr) cos 0 + A GAr  + dr)sinG  + Az. 

rx-r2  y 


We  need  to  find  the  Max  and  Min  of  / (r,  6)  - -Jx' 


2+/2 


R — R 

-J z , to  eliminate  r and  6.  However, 


we  cannot  yet  solve  the  two  non-linear  equations  formed  by  equating  the  partial  derivatives  of  the 
function  / with  respect  to  r and  t to  zero  in  compact  (symbolic)  form.  So  in  order  to  use  the  results 
for  the  conical  surfaces  in  the  tolerance  synthesis/assemblability  analysis,  we  have  to  explore  further  to 
get  a solution  or  else  develop/adopt  methods  to  treat  the  equation  keeping  the  two  free  parameters  as 
extra  variables. 


Page  26 


Tolerance  Synthesis  Scheme 


Page  27  of  60 


3.3.7  Special  cases 
3.3.7.1  Composite  tolerance  specification 

Sometimes  we  have  more  than  one  tolerance  specified  on  one  feature  of  size,  as  shown  in  Figure  19.  In 
this  case,  more  constraints  need  to  be  generated.  A designer  has  two  intents  when  specifying  combined 
tolerances:  1)  both  tolerances  must  be  satisfied  at  the  same  time,  and  2)  some  DOFs  of  the  feature  of 
size  will  be  constrained  specifically  by  the  additional  tolerance.  In  the  following  case,  the 
perpendicular  tolerance  gives  more  constraint  on  the  orientation  of  the  feature  of  size. 
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Figure  19:  Composite  tolerance  specification 


The  ASME  Y14.5  standard  requires  that  all  the  tolerances  specified  to  one  feature  of  size  must  be 
satisfied  independently.  This  means  that  we  need  to  treat  every  tolerance  seprately.  For  the  above 
example,  the  constraints  for  positional  tolerance  are  the  same  as  those  in  Section  3.3.3;  the  following  is 
the  derivation  of  the  constraints  for  perpendicular  tolerance. 

The  function  of  the  perpendicular  tolerance  here  is  to  give  more  constraints  on  the  orientation  of  the 
feature.  Because  perpendicular  tolerance  does  not  fix  the  position  of  virtual  condition  boundary,  and 
the  tolerance  only  controls  orientation,  the  position  of  the  transformed  point  due  to  the  perpendicular 
tolerance  can  be  calculated  as  follows: 


1 

0 

My 

o' 

(r  + dr)  cosd 

(r  + dr)  cos#  + A 0vz 

0 

1 

-a  ex 

0 

(r  + dr)  sin  ^ 

(r  + dr)  sin  6 - A0xz 

-a  e 

y 

<3 

i 

0 

z 

z - A0y(r  + dr) cos#  + A0x(r  + dr)  sin 0 

0 

0 

0 

1 

1 

1 

Therefore,  the  constraint  equation  will  be: 

((r  + dr)  cos  d + A#vz)2  + ((r  + dr)  sin  <9-  A0xz)2  < ( 2r  + TU  + TV)2  / 4, 
with  the  constraint  equation  by  positional  tolerance: 

((r  + dr)  cos  d + A0 z + Ax)2  + ((r  + dr) sin#  - A 0xz  + Ay)2  < (2r  + TU  + Tp)2  /4 


(3.3.7. 1-1) 


(3.3.7. 1-2) 
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Equation  (3.3.7. 1-1)  and  (3.3.7. 1-2)  are  the  constraints  due  to  the  composite  tolerance  specification  in 
the  above  example.  Using  the  same  procedure  as  in  section  3.3.3  to  eliminate  6 and  z,  we  have 
thefollowing  constraints  on  deviations  for  the  perpendicular  tolerance: 
f'(9,z ) = ( a cos  9 + bz)2  + (a  sin  9 + dz)2 , 
where  a = r + dr,  b = A9v,  d = -A9X. 


By  equating  the  partial  derivative  s of  f ' to  zero, 

= 0,  and  - 0,  we  get, 

89  dz 


fL  = max 


±1 


(ad2  + ab2) 


d2  +b- 


/ ' = min 

J min 


+i 


Vrf2+A: 


(ad2  + ab2) 


d2  + bA 


From  the  last  section  for  a positional  tolerance  of  similar  specification,  we 

have: 

/ (6,z)  = (a  cos  9 + bz  + c)2  + (a  sin  0 + dz  + e)2 , 

where  a = r + dr,  b = A9y  , c-  Ax,  d--A6x,  e = Ay. 

By  equating  the  partial  derivative  s of  f to  zero, 

-f—  = 0,  and  — = 0,  we  get, 

89  dz 


/max  = max 


±1 


yjd2  +b: 


(ad2  + ab2)  - be  + cdy 


d2  +b< 


-,o 


»,  = min 

5 J min 


±1 


\Jd 2 +b: 


:(ad2  + ab2)  - be  + cd)' 


d2  +b J 


-,0 


The  final  constraints  due  to  perpendicular  tolerance  and  positional  tolerance  are: 


f < 

J max 


f < 

J max 


(2r  + TU  + Tp 


2 r + TU  + Tv 


(3.3.7.1-3) 


-TL<  dr  < TU . 
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When  the  toleranced  feature  is  internal,  as  for  example,  for  a hole,  the  constraints  are: 

r2r  + TU-Tp 


f > 

J min 


f ^ 

J min  — 


V 


(2r  + TU -Tv 


-TL<  dr<TU. 


(3.3.7. 1-4) 


If  the  external  feature  is  modified  to  LMC,  the  constraints  will  be: 

f 2r-TL-Tp^ 


f > 

J min 


f ^ 

J min 


v 


2r  -TL-  Tv 


-TL<  dr  < TU . 


If  the  internal  feature  is  modified  to  LMC,  the  constraints  will  be: 


f < 

J max 


f ^ 

J max  — 


2r-TL  + TpV 

2 J ’ 

f 2 r-TL  + TvV 


V 2 j 

-TL<  dr  < TU . 


3.3.7.2  MMC  Specification  on  Datums 


<f>  2 r±% 


(3.3.7. 1-5) 


(3.3.7. 1-6) 


Figure  20:  MMC  Specified  to  datum 

When  accurate  relative  position  of  two  parts  is  not  required,  the  designer  may  specify  MMC  on  the 
datum  to  allow  more  relative  position  variation  between  a datum  feature  of  size  and  the  toleranced 
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feature  of  size,  if  the  size  of  the  datum  feature  of  size  departs  from  its  MMC.  This  specification  will 
not  increase  the  size  of  the  VCB  but  will  allow  its  center  axis  to  move  if  the  datum  feature  of  size 
departs  from  its  MMC. 

In  contrast  to  the  composite  tolerance  specification  described  in  the  previous  section,  the  position  of 
the  VCB  of  the  left  hand  side  cylindrical  feature  is  dependent  on  the  real  size  of  the  cylindrical  feature 
on  the  right  hand  side.  Because,  this  specification  means  that  both  cylindrical  features  will  mate 
features  on  the  other  part.  The  clearance  due  to  the  departure  of  size  of  datum  feature  of  size  from  its 
MMC  is  allowed  to  be  used  to  adjust  the  position  of  left  side  cylindrical  feature,  this  problem  can  be 
modeled  as  if  we  have  the  positional  tolerance  and  perpendicular  tolerance  specified  on  the  left-hand 
feature. 

The  positional  tolerance  will  be  Tpi+(TU+Tp-dr)  6 and  the  perpendicular  tolerance  is  Tp].  Following 
the  procedures  as  described  in  Section  3.2,  we  can  develop  the  constraints  on  deviations  of  the  left 
hand  side  cylindrical  feature  by  substituting  Tp=  Tpi+(TU+Tp-dr)  and  Tv=Tpi : 


f < 

J max 


f ^ 

J max 


2r  + TU  + Tpl  + (TU  + Tp  - dr) 


2r  + TU  + Tp\ 


-TL<  dr  < TU . 


Example  of  Some  Mappings 
Case  1:  Planar  Feature 


Figure  21:  Planar  feature  (rectangular)  with  size  tolerance 


6 In  this  case  ( TU+Tp-dr j is  from  the  other  cylindrical  feature 
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Figure  21  shows  a simple  size  tolerance  to  a planar  surface  with  LCS  specified.  Let  us  assume  the  left 
hand  side  surface  is  the  implicit  datum  and  the  other  side  is  controlled  by  the  tolerance  specification. 
From  Figure  21,  we  know  the  deviation  parameters  for  the  right  hand  side  surface  which  is  a non-size 
feature  are  A0X  A6y,  Az. 

The  following  constraints  are  established  by  considering  the  deviation  at  the  four  extreme  points 
(comers)  of  the  plane  which  are  nominally  at  (1.5,  2,  0),  (1.5,  -2,  0),  (-1.5,  2,  0)  and  (-1.5,  -2,  0). 


-0.08  < Az  + 1.5*A0y- 2*A6X  < 0.08, 

-0.08  < Az  + 1.5*A6y  + 2*A6X  < 0.08, 

-0.08  < Az-  l.5*A0y  + 2*A6X  < 0.08, 

-0.08  < Az-  1.5*A6y  - 2*A0X  < 0.08. 

The  above  constraints  on  deviation  parameters  result  in  a diamond  shape  in  deviation  space  as 
follows:. 


Figure  22:  Deviation  space  for  size  tolerance  on  planar  surface 


Case  2;  Cylindrical  mating 


Figure  23:  Assembly  of  two  parts  specified  with  positional  tolerance 
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From  the  above  specification,  we  know  that  the  VCB  of  the  external  feature  of  the  above  part  is  equal 
to  that  of  the  internal  feature  of  the  bottom  part.  Following  equations  (3. 3. 3-1)  and  (3. 3. 3-2),  the  size  of 
the  VCB  for  the  external  cylindrical  feature  is  5+0.25+0.05=  5.3,  and  that  of  the  internal  cylindrical 
feature  is  .5.7-0.35-0.05=  5.3.  Therefore,  we  know  that  the  two  parts  are  guaranteed  to  be  assembled. 
Now  let’s  see  what  constraints  on  deviation  parameters  we  can  get  by  our  tolerance  mapping,  so  we 
can  computationally  perform  the  tolerance  analysis  or  synthesis. 


Let  a = r + dr,  b - AO  , c = Ax,  d = -A0x,  e = A y, 


= max 

•J  max 


±1 


yjd2  +b: 


:(ad2  + ab2)  - be  + cd 


d2  +b‘ 


,0 


, = min 

’ J min 


±1 


V 1 


Jd2+b: 


(ad2  + ab2)-  be  + cd 


d2  +bl 


-,o 


The  final  set  of  constraints  for  the  external  feature  is: 

, ( 5 + 0.25  + 0.05^2 

f < 

J max 


V 


2 j 

-0.25<  dr  < 0.25. 


When  the  toleranced  feature  is  internal,  as  for  example,  for  a hole,  the  constraints  are: 

f > 

J min  — 


^5.7  -0.35  -0.05^2 


-0.35<  dr  <0.35. 

The  above  equations  cannot  be  reduced  further  in  compact  form  so  no  further  derivation  and/or 
plotting  are  possible.  However,  the  inequalities  as  such  could  be  used  in  computations. 


4.  Optimization  Process 


Formulation  of  the  Cost  Functions 

The  cost  function  to  be  minimized  can  have  several  components,  including: 

o Cost  of  Manufacturing, 
o Cost  of  Assembling, 
o Cost  of  Inspection, 
o Cost  related  to  product  life  span, 
o Cost  to  meet  functional  requirements, 
o Cost  to  meet  assemblability,  etc. 

In  this  study,  we  have  considered  the  Cost  of  Manufacturing  in  details  and  a new  cost  of  manufacturing 
model  [deviation-based  cost  of  manufacturing]  has  been  introduced.  Apart  from  the  cost  of 


Page  33 


Tolerance  Synthesis  Scheme 


Page  34  of  60 


manufacturing,  we  could  not  yet  carry  out  further  investigation  into  the  cost  of  assembly,  inspection, 
product  life-span  related  costs,  etc.  Some  authors  have  done  context  specific  work  on  the  above  costs 
and,  for  example,  Dong  [DON97]  has  considered  the  effect  of  tolerance  on  cost  related  to  product 
lifespan  by  constructing  a functional  performance  index  generated  from  the  effect  of  tolerance. 
However,  the  work  is  applicable  to  very  specific  domains  and  cannot  be  used  universally. 

In  our  tolerance  synthesis  scheme,  we  have  not  considered  these  costs.  In  the  proposed  tolerance 
synthesis  scheme,  we  have  kept  a generalized  cost  of  manufacturing  model  so  that  the  designer  can 
easily  mtroduce/modify  the  cost  functions  associated  with  each  part  in  the  artifact  library  [ROYOl]. 
The  cost  function  has  been  kept  as  a virtual  function  in  the  artifact  class  so  that  the  specific  cost 
models  could  be  introduced  in  the  instantiated  object  from  a class  of  real  cost  functions. 

We  would  consider  the  last  two  of  the  above  cost  factors  (Cost  to  meet  functional  requirements  and 
Cost  to  meet  assembleability)  as  constraints  rather  than  minimization  criteria.  We  will  now  go  into 
details  of  cost  of  manufacturing. 


4. 1 Cost  of  Manufacturing 

We  consider  the  cost  associated  with  manufacturing  a part  for  specified  material, 
dimensions/geometric  shape,  a sequence  of  manufacturing  operations,  and  the  process  capabilities.  In 
general,  the  cost  of  manufacturing  is  a function  of  all  the  above  entities.  Many  researchers  have 
analyzed  various  issues  associated  with  the  cost  of  manufacturing  in  relation  to  tolerances  and  the  cost 
of  manufacturing  as  a decreasing  function  of  a single  tolerance  parameter.  Various  functional  forms 
(like  the  inverse  power  law,  exponential  decay,  etc.)  have  been  used  as  cost  models  that  are  generally 
monotonically  decreasing  functions  of  increasing  in  tolerance  value.  Effect  of  process  capabilities  on 
cost  of  manufacturing  has  also  been  well  studied  by  several  authors  [DON97], 

Also,  since,  in  general,  more  than  one  operation  is  required  to  transform  the  raw  blank  into  the  final 
finished  part,  the  cost  of  manufacturing  is  a function  of  the  process  sequence  and  of  how  much 
accuracy  is  achieved  in  each  stage  of  operation.  The  cost  is  also  affected  by  the  setup  error  in  each 
machining  process.  The  total  cost  of  production  thus  becomes  a sum  of  the  costs  associated  with  each 
process  [ROY97]. 

While  none  of  these  methods  mentioned  above  could  claim  to  be  universal,  there  are  several 
limitations  with  the  one-parameter  (single  tolerance)  cost  of  manufacturing  formulation.  In  reality,  a 
manufactured  surface  would  rarely  have  a single  tolerance  value.  Apart  from  a size  tolerance,  there 
could  be  geometric  tolerances  (form  tolerance,  orientation  tolerance,  etc.)  and  it  would  be  difficult  to 
formulate  a single  parameter  representing  all  these  tolerances  that  could  effectively  be  used  for 
representing  the  cost  of  manufacturing.  In  order  to  circumvent  the  above  problem,  we  propose  a new 
formulation  for  the  cost  of  manufacturing,  which  we  would  like  to  call  the  deviation-based  model. 


4.2  Deviation-Based  Cost  Model  of  Manufacturing 

In  general,  a part  will  have  more  than  one  feature,  and  is  connected  to  other  features  of  other  parts  and 
since  each  feature  could  possibly  have  up  to  six  degrees  of  freedom,  the  scheme  we  have  adopted  uses 
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six  deviation  parameters  to  represent  the  variations  associated  with  each  feature.  In  our  earlier  reports, 
we  have  used  the  small  deviation  torsors  (SDT)  for  representing  the  deviations  of  features.  We  have 
also  shown  that  these  deviation  parameters  are  related  to  the  tolerance  specifications  through  a 
functional  mapping  (a  one  to  many  relation:  a tolerance  specification  defined  a bounded  region  in  the 
deviation  space).  In  our  approach,  we  propose  to  represent  the  cost  of  manufacturing  as  function  of 
these  deviation  parameters.  Thus,  for  example,  the  cost  function  g(S)  defined  as  a function  of  some 
tolerance  value  8,  would  become  a function  of  the  six  parameters  g(S)  = g(  8(0X  9V  9Z  5x  5y  5Z),  and 
then  the  cost  functions  could  be  treated  as  functions  of  the  small  deviation  parameters  associated  with 
the  surface/part. 

We  propose  to  model  the  cost  of  manufacturing  a part  as  an  explicit  product  of  six  functions  of  the  six 
deviation  parameters  in  the  form: 

DCOM(d)  = Ci(d1)*C2(d2)*....C6(d6),  (4.2-1) 

where  d = (di,  d2,  d3,  d4,  ds,  df,)  = (9X,  9y,  9- , 8X , 8,, , 8J  is  a 6-component  torsor  representing  the  six 
deviation  parameters  characteristic  of  the  feature. 

The  form  of  the  individual  functions  can  vary  depending  on  the  specific  surfaces  and  manufacturing 
process.  Some  of  the  functions  could  be  of  similar  form.  Also  depending  on  the  nature/type  of  the 
feature,  some  of  the  functions  will  be  constants  (i.e.,  Vx,  f(x)  = 1 ) and  could  be  eliminated.  These 
functions  correspond  to  the  deviation  parameters  that  are  invariants  of  the  surface. 

For  an  example,  for  a planar  surface  there  are  only  three  independent  parameters  given  by:  d = (dx,  0, 
0,  0,  0V,  9Z)  that  affect  the  deviation  of  the  surface  from  it’s  nominal  shape;  changes  in  the  remaining 
three  keep  the  surface  invariant.  The  cost  function  can  then  be  represented  as: 

DCOM(d)  = Cx(dx)  *Cy 0y)  *Ce(0z).  (4.2-2) 

Also,  in  this  case,  since  the  surface  is  symmetric  about  the  y and  z axes,  the  form  the  two  functions  for 
rotational  deviations  along  these  two  directions  are  also  the  same,  namely  Cq. 

We  have  mentioned  earlier  in  this  section  that  the  deviation  parameters  are  restricted  by  the  tolerance 
specification.  For  the  planar  feature,  this  mapping  forms  a convex  hull  in  the  form  of  a diamond  in  the 
3-D  space.  For  a rectangular  planar  section  with  cross-section  ( ‘a’  x ‘b’),  (see  section  3.3.2),  the 
following  are  the  constraints: 


Tsl  ^ min  (dx+  a*0y  + b*  0Z  , dx+  a*0y  - b*  6Z , dx-  a*9y  + b*  0Z , dx-  a*6y  - b*  6Z,  (4.2-3) 

Tsu^max  (dx+  a*9y  + b*  6Z , dx+  a*0y  - b*  9Z  , dx-  a*0y  + b*  9Z,  dx-  a*0y  - b*  9Z),  (4.2-4) 

where  (Tsl,  Tsu)  are  the  lower  and  upper  values  of  the  tolerance  parameter  for  the  planar  surface. 
Thus,  the  parameters  dx  9y  9Z  are  restricted.  The  cost  function  (3-2)  will  be  restricted  accordingly. 

To  illustrate  the  cost  function,  let  us  assume  a generic  function  of  the  form:  C(x)  = a +b/ |x|,  where  x is 
the  deviation  parameter  and  a and  b are  constants.  For  the  planar  case,  we  then  get, 

DCOM(d)  = Cx(dx)*Cg(9y)*Cg(0z)=(ai+bi/\dxO  * (a2+b2/\  9y()  * (a3+b3/\0z\) 
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To  further  simplify,  we  assume,  a=b=l  and  a 2 = 03,  b2  = by 
Then  we  have 


-Tsl  ^ (dx+  0y  + 0z)  ^ Tsu ; 


- Tsl  ^ (dx+  %>  - &)  ^ Tsu . 

-Tsl  ^ (dx-  Gy  + 0Z)  <Tsu, 

- Tsl  ^ (dx-  Gy  - Gz)  -S' Tsu . 

and  DCOM(d)  = (a,+b,/\dx\)  * (a+b/\  Gy\)  * (a+b/ 1 0Z\). 

Removing  the  z-parameter  6Z  again  so  that  we  can  have  a visual  representation  of  the  cost  function,  we 
have  - Tsl  -S'  (d+  G)<Tsu  -Tsl  - (d-  6)  <Tsu  and 

DCOM(d)  = (a,+b,/\d\)  * (a+b/\  6\).  (4.2-5) 

In  the  d-6  plane,  this  would  look  like  a tent  bounded  by  four  vertical  planes  (by  the  tolerance 
specification)  approaching  infinity  along  the  two  axes  (Figure  24,  Figure  25). 


Cost 


Figure  24:  Cost  as  function  of  deviation  parameters 
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- TSl  ^ (d-  0)  <Tsu 


e 


-Tsl  < (d+  6)  <Tsu 


Figure  25:  Cost  contour  lines  and  the  bounds 


It  is  to  be  noted  that  the  actual  form  of  the  individual  cost  functions  are  yet  to  be  determined.  Also,  it  is 
evident  that  since  the  deviation  parameters  limit  the  tolerance  values  in  some  way  or  other,  the  cost 
functions  of  these  deviation  parameters  should  also  be  monotone  non-increasing  in  nature. 


4.3  Optimization  Process 

Based  on  the  models  described  above,  the  optimization  process  could  be  summarized  as  below: 
Minimize: 

Cost  of  production.  C(x)  — F(j  e all parts)E(je  all  mating  features)  F(Xjj),  (4.3-1) 

where  j is  the  deviation  parameters  of  the  j feature  of  part  / in  the  assembly 
Subject  to: 

Functional  constraints  of  the  form:  G( x)  <0,x  e R"  , (4.3-2) 

Linear  Assembleability  Constraints  of  the  form:  Ax  = b , where  A = /a,y/  is  an  m x n matrix, 
dij€  R,  x e Rn  b e Rm  , and  m is  the  number  of  constraints.  In  general,  m is  less  than  the 
number  of  variables  n. 

Deviation  to  tolerance  mapping  constraints:  FDPk(x)  < FTolk(t),  xeRn,  ke  N*  (4.3-3) 
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where  FDP  denotes  Function  of  Deviation  Parameter,  FToI  denotes  Function  of  tolerance  parameter 
subscript  k is  the  number  of  constraint  equations. 

In  order  to  look  for  a solution  for  the  above  minimization  formulation  that  falls  under  the  non-linear 
programming  [NLP]  category,  we  are  proposing  to  develop  suitable  and  efficient  evolutionary 
algorithms  [EA]  because  of  the  potential  problems  that  could  be  anticipated  in  using  standard  gradient- 
based  NLP  solvers  as  detailed  below. 

The  problem  formulated  above  has  several  layers  of  complexity  that  make  the  problem  very  difficult  to 
solve  using  traditional  methods.  We  have  tried  to  formulate  the  constraints  for  a very  simple, 
elementary,  two  dimensional,  planar,  3-element  assembly  [Appendix-A]  for  which  we  need  about 
150+  variables,  80+  linear  constraints,  10+  non-linear  constraints.  We  assumed  there  were  no 
constraints  arising  out  of  the  function  requirements  and  assumed  a fairly  complex  [exponential  and 
inverse]  cost  function.  We  then  tried  to  solve  the  problem  using  the  Matlab©  NLP  solver  fmincon.  We 
got  some  solutions  after  several  trials  but  the  solution  is  not  realistic  unless  we  impose  several 
restrictions  on  the  deviation  parameters  as  well  as  the  tolerance  parameters.  We  anticipate  that  a 
medium  to  reasonably  large-sized  assembly  with  say  a couple  of  hundred  parts  in  a 3D  space  would 
give  rise  to  several  thousands  of  variables  and  constraints.  Apart  from  that  the  matrix  A in  the  linear 
assembleability  constraints  Ax  = b could  be  very  sparse  depending  upon  the  assembly,  requiring 
special  treatment. 

In  addition,  it  is  our  endeavor  to  establish  a deviation  hull  as  a suitable  solution  space  for  further 
mapping  of  the  tolerance  space.  This  task  itself  is  quite  complicated  and  we  believe  an  approach  using 
the  EA  techniques  could  give  us  some  lead  by,  for  example,  storing  the  intermediate  solutions  and 
doing  some  analysis  of  the  pattern.  However,  this  area  is  open  for  further  research. 


Computer  Implementation  of  the  Tolerance  Synthesis  System 

The  proposed  tolerance  synthesis  system  has  been  considered  as  an  object-oriented  representation  built 
in  part  on  our  earlier  work  [ROYOl]  where  we  introduced  a generic  artifact  class  that  can  be  extended 
to  include  a tolerance  class  and  other  class  structures.  We  developed  the  class  structures  using  UML 
(Rational  Rose  2000©  - Enterprise  Edition).  We  plan  to  do  the  full-scale  implementation  in  C++  and 
add  interfaces  with  3D  CAD  packages  such  as  AutoCAD©,  ProE©,  and  SolidWorks©  etc.,  so  that  the 
end  user  could  interact  with  a user-friendly  graphic  environment  for  defining  the  product  and  carry  out 
the  desired  tolerance-related  computations.  At  present  the  graphics  display  and  other  graphics  driven 
design/tolerance  synthesis  aspects  are  being  formulated  for  use  with  SolidWorks©.  Appendix  - B 
shows  a part  of  the  class  diagrams  of  the  basic  entities  in  the  synthesis  scheme. 

A very  skeletal  (but  operational)  tolerance  synthesis  method  has  been  implemented  under  the 
Windows©  environment  as  a core  module  without  Windows  GUI  and  CAD  interfaces  using  a database 
programming  system  (Clipper  5.2©).  In  order  to  make  manipulation  of  the  model  data  easy, 
systematic,  and  structured,  the  entire  data  associated  with  each  model  is  stored  in  a relational  database 
system  with  several  tables  for  each  type  of  object/class.  The  structures  of  these  tables,  along  with 
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relational  links  between  the  tables,  have  been  shown  in  Appendix-C.  The  field  names  in  each  table 
have  been  kept  nearer  to  their  real  counterparts  for  ease  of  understanding. 

The  program  is  menu  driven  and  operates  in  two  phases.  In  the  first  phase,  the  designer  would  build  the 
model  (artifact)  using  standard  artifact  library  elements  and  specifying  other  details.  After  a model  has 
been  established,  governing  equations  for  each  independent  path  are  generated  and  then  the  equations 
are  converted  into  a set  of  MatLab©  input  files  for  optimization.  Cost  functions  and  constraints  for  the 
functional  requirements  are  formulated  separately  as  MatLab©  input  files.  MatLab©  is  then  used  to 
solve  the  resulting  nonlinear  optimization  problem  for  arriving  at  a set  of  optimal  deviation  parameters. 

Details  of  the  Example  Model 

In  Appendix  - D we  present  details  of  an  example  model  including  part/feature  details,  input  data 
listing,  generated  data  (partial  list),  generated  MatLab©  input  files  (mostly  generated  from  the 
synthesis  module),  and  some  preliminary  results.  It  has  already  been  mentioned  that  we  are  using 
MatLab©  to  solve  the  cost  minimization  problem.  In  the  present  case,  only  a sample  (fictitious)  cost 
function  has  been  used  to  test  if  the  constraints  generated  could  be  properly  treated  in  MatLab©.  The 
optimized  output  results  are  of  little  significance  here  since  those  are  based  on  fictitious  cost  functions. 
We  are  now  in  the  process  of  selecting  realistic  cost  functions  and  various  models,  so  that  some 
practical  tolerance  limits/bounds  could  be  computed  based  on  the  optimization  results. 
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Appendix  - A:  An  Example  of  Deriving  Assembleability 

Constraints 


Mathematical  Model  of  Tolerances 

Various  mathematical  models  for  representation  of  components  in  an  assembly  have  been  studied  so 
that  the  geometric  variations  arising  out  of  small  deviations  of  the  features  (faces)  of  the  components 
could  be  analyzed.  We  adopted  the  small  displacement  torsor  scheme  to  represent  the  deviation  of  a 
real  surface  with  respect  to  the  corresponding  nominal  surface  at  a point  on  a feature/face  of  a 
component.  Based  on  this  scheme,  an  effective  representation  of  kinematic  variations  for  the  mating 
faces/features  could  be  formulated  as  elaborated  hereinafter. 

A torsor  is  a tensor-type  representation  of  small  displacements  and  rotations  (corresponding  to  the  six 
degrees  of  freedom  associated  with  a space  point  on  a surface  embedded  in  R3).  If  we  introduce  three 
small  displacements  along  the  three  local  principal  axes  at  a point  A on  a surface  as  8 = [S*  5y  SZ]T 
and  three  rotations  0 = [6X  0y  6Z  ]T  then  the  torsor  Ta/r  at  A with  reference  to  a datum  R can  be 
represented  as: 


T'a/r  ~ [9a  §a  ]T  ~~  [ Ox  &y  Oz  8X  Sy , Sz  ]T. 


Torsors  transform  according  to  standard  rotation  and  translation  rules  of  axes.  However,  the 
components  are  small  in  comparison  to  the  dimension  of  the  components. 

The  torsor  representation  of  deviations  of  real  surfaces  with  respect  to  the  nominal  surface  for  all 
possible  regular  surface  geometries  used  in  engineering  design  have  been  well  studied.  It  has  been 
found  that,  there  can  be  some  component  (out  of  the  possible  six)  variations  along  which  the  surface  is 
kept  invariant  with  respect  to  the  nominal  surface.  These  components  are  called  independent 
components  and  may  be  varied  for  considering  possible  configurations.  For  example,  for  a cylindrical 
surface,  a translation  along  the  axis  keeps  the  surface  invariant  as  does  a rotation  about  the  axis.  Thus 
the  torsor  representation  for  the  deviation  is: 


Ta/r  ~ [&a  &a]  — [Rx  0y  6Z  lx  8y  Sz] 
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where  / stands  for  independent  movement  (in  this  case  along  the  axis  of  the  cylinder)  and  R stands  for 
free  rotation  (here,  about  the  axis  of  the  cylinder).  We  can  arbitrarily  make  these  quantities  zero  so 
that  they  do  not  enter  into  the  equations.  In  addition,  we  can  vary  these  parameters  for  various 
constraint  satisfactions  keeping  the  surfaces  invariant. 

The  small-displacement  torsor  mentioned  above  could  be  used  to  represent  three  types  of  deviations 
associated  with  the  mating  of  components  in  an  assembly.  Two  torsors  associated  with  each 
component  and  a gap  torsor  at  the  mating.  The  two  torsors  associated  with  each  surface  are:  a 
displacement  torsor  and  a deviation  torsor  as  shown  in  Figure  26.  The  displacement  torsor  is  an 
intrinsic  variation  torsor  that  is  specific  to  the  element  (for  example,  displacement  arising  out  of  the 
variation  of  the  radius  of  a cylindrical  surface).  The  displacement  torsor  represents  the  deviation  of  the 
real  surface  from  the  nominal  surface.  The  gap  torsor  between  two  mating  features  represents  possible 
geometric  variations  between  the  mating  features,  as  depicted  in  Figure  26. 


Figure  26:  Interacting  torsors  at  mating  surfaces 


The  torsor  scheme  indicated  above  could  be  used  to  form  equations  based  on  loops  connecting 
components  in  an  assembly.  Torsors  associated  with  each  component  feature  would  be  kept  in  the 
artifact  library  in  the  local  coordinate  system  of  the  component.  The  gap  torsor  associated  with  two 
mating  features  would  be  formed  based  on  the  functional  and/or  kinematic  requirements  as  established 
in  the  design  specifications.  In  order  to  add  the  torsors  in  a loop  to  form  equations  connecting  the 
components  of  the  torsors,  a mechanism  would  be  necessary  to  transform  the  torsors  defined  in  the 
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local  coordinate  systems  of  the  individual  components  into  a global  reference  frame.  Since  a torsor  is 
represented  by  a 1x6  vector,  the  transformations  required  would  be  6x6  transformation  matrices. 

Torsor  Transformation  Rules 

T T 

The  small  displacement  torsor  Tn  = [On  $n]  = [0X  Qy  0Z  & Sy  &]  defined  at  a point  iVona  surface 
could  be  expressed  at  a point  M on  a different  surface  by  applying  the  required  displacement  and 
rotation  to  the  vector  MN.  The  transformation  is  a 6x6  matrix  given  by: 


1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

- Az 

Ay 

1 

0 

0 

Az 

0 

-Ax 

0 

1 

0 

-Ay 

Ax 

0 

0 

0 

1 

where  (Ax  Ay  Az)T  is  the  vector  from  M to  N. 

The  transformation  is  given  by  Tm  = DTn  and  it  can  be  explicitly  written  as 

0 x 
0 y 

T - 

M S x + (A y6  z - A z6  y 
8 y + (tx  z 0 x - Ax  0 z 
8 z + (A  x 0 y - A y 0 x 


0 


N 


8 N + MN  x 0 N 


The  torsor  scheme  for  the  modeling  of  the  small  geometric  variations  of  connecting  features/faces  of 
an  artifact  has  been  elaborated  with  a simple  example  of  a rectangular  block  (Part-3)  enclosed  within 
two  channels  (Part-1  and  Part-2)  (Figure  27,  Figure  28).  Apart  from  the  simplified  configuration, 
further  simplification  has  been  done  as  detailed  below  so  that  the  relevant  transformations  and 
preliminary  computations  could  be  checked  manually. 
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Figure  27:  Example  problem  with  three  blocks 
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Figure  28:  Example  problem  with  feature  numbers  and  LCS 

The  simplifications  are: 

a)  All  variations  take  place  along  the  vertical  plane  only  (with  3 possible  variations:  x- 

rotation,  ^-displacement  and  z-displacement)  to  keep  the  formulation  simple. 
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b)  There  is  no  displacement  torsor  for  each  part.  Only  the  surface  deviation  torsors  and 
the  gaps  are  taken  into  consideration. 

c)  The  bottom  surface  of  Part  1 has  been  assumed  to  be  a fixed  datum. 

The  three  parts  1,  2,  and  3,  local  axes  and  contact  surfaces  have  been  shown  in  the  same  sketch. 
Surfaces  have  been  designated  as  <part  number><seria!  number>  such  that  surfaces  on  part  1 are 
numbered  as  11,  12,  13,  14,  etc.  For  generating  equations  connecting  the  variations  of  the  mating 
surfaces,  six  loops  have  been  considered  (Figure  27). 

With  the  above  simplifications,  we  are  left  with  following  configuration: 


Number  of  mating  surfaces 
Number  of  gaps  between  mating  surfaces 
Number  of  surface  parameters 
Number  of  gap  parameters 
Total  number  of  parameters 
Number  of  loops 

Number  of  equations  from  the  loops 

The  six  loops  are  as  follows: 

Loop-1:  12-^32-^31-^13 
Loop-2:  1 1->21->22->32->32->  12 
Loop-3:  22-^32-^33-^23 
Loop-4:  24-^34->33~>23 
Loop-5:  15->25-^24-^34->34->  14 
Loop-6:  14-^34->31-^  13 

The  associated  torsors  in  their  local  coordinate  systems  are: 


m tan  pianer; 
8 

14x3  =42 
8x3  =24 

42  + 24  = 66 
6 

6x3  = 18 


rw~'  rri  rr  rr  rr  r t ' ri ' r 1 ' rr  rtr  rwy 

illy  i 12y  i 13y  i 14y  i 15j  i 21y  * 22y  * 23y  * 24y  * 25j  * Sly  * 32y  * 33y  *34. 

Typically,  each  of  these  torsors  would  have  similar  6-component  forms  like: 
Tn  =[0u  Su  [T  ~ [0xii  0 0 0 Syu  SzuJT  , 

T 12  -[012  012  f - [0x12  0 0 0 Syi2  0zl2  [T‘ 


The  gap  torsors  are: 

G 11/21,  G 12/32,  G 13/31,  G 14/34,  G 15/25,  G 22/32,  G 23/33,  G 24/34, 
with  typical  forms  given  by: 

T T 

G n/21  = [0g  11/21  0§  11/21  ] =[0gxll/21  0 0 0 Sgyii/21  dgzu/2l  J . 
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We  now  consider  the  Loop  - 1 (12->32->3 1 ->  13)  and  write  the  torsor  equation  for  the  loop  as: 


T]2  + G 12/32  + (-T32)  + T31  + (-G 13/31 ) + (-T13)  - 0. 


Transforming  all  the  terms  in  the  above  equation  in  Part-l’s  frame  of  reference,  the  equation  becomes, 

T]2  + (D 12/23)  12  G 12/32  + (1)32)12  (-T32)  +(T>3l)l2  T 31  +(T>l3/3l)l2  G 13/31  + (D  13)12  (-T13)  = 0, 

where  the  D ’s  are  the  corresponding  transformations  with  respect  to  the  Part  - 1 . 

Writing  the  D ’s  in  full,  and  carrying  out  the  operations,  the  above  equation  becomes, 

10X12  0 0 0 Syl2  8*2  f + 

/ Ogx 12/32  0 0 0 8gyi2/32  Sgz  12/32  / ~ 

[6x32  0 0 0 Sy32  + 6x32  (Z32~  Z12)  &J2  + 0x32  (y32~  J 12)  f + 

/ 0x3i  0 0 0 Sy3i  + 6x31  (Z31  - Z12)  8*1  + 6X3i  (y>3i  - yi2)  1T  + 

[0gxl3/31  0 0 0 Sgy  13/31+  6gx  13/31  (Z31~  Z12)  Sgzl3/31  + 0gxl3/31  (y>31  ~ y 12)]T  “ 

[6x13  0 0 0 Syi3  + 6x13  (Zl3-  Z12)  8zl3+  0xl3(yi2- yi2)]T 

= [0  0 0 0 0 0]T. 


Equating  each  component  to  zero  gives  six  scalar  equations.  In  this  present  case,  there  are  only  three 
equations  due  to  the  restricted  space  of  the  example. 

6x12  0gxl2/32  ~ 0x32  0x31  + 0gxl3/31  ~ 0x13  0, 

Syi2  + Sgy 12K2  ~ Sy32  - 0x32  (Z32  ~ Z12)  + Sy31  + 0x31  (Z31  ~ Z12)  + Sgy  13/31  + 

0gx  13/31  (Z31  - Z12)  ~ Syi  3 + 0x13  (Zl3  ~ Z12)  ~ 0 , 

SZ12  + Sgzl2/32  - Sz32  ~ 0x32  (y32  ~ yi2)  + S*i  + 0X31  (y31  ~ yi2)  + Sgz]3/3i  + 

Ogx  13/31  (y3i~yi2>  ■ Sri 3 + 0xi3  (yn-yi2)  = o. 


Each  of  the  six  loops  will  give  three  similar  equations.  Altogether  eighteen  equations  would  be  formed 
connecting  the  sixty-six  parameters. 

Let  us  repeat  the  procedure  of  forming  the  equations  for  Loop-2  (1 1->21->22->32->32->  12). 

The  torsor  equation  for  the  Loop  - 2 is: 

T11  + G u/21  + (-T21)  + T22  + (G 22/32)  + (-T32)  + T32  + (G 32/1 2)  + (-T12)  — 0, 
i.e.,  T11  + G n/21  + (-T21)  + T22+  (G 22/32)  (G32/12)  + (-T12)  = 0. 
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Transforming  all  the  terms  in  the  above  equation  in  Part-1  ’s  frame  of  reference,  the  equation  becomes, 

Tu  + (D ll/2l)ll  G 11/21  + (D 2 l)ll  (-T21)  +(D22)ll  T22+(D 22/32)  11  G 22/32 
+ (D  12/12)11  G22/12  + (Did  11  (-T12)  = 0, 


where  the  D ’s  are  the  corresponding  transformations  with  respect  to  the  Part- 1 . 
Writing  the  D’s  in  full,  and  carrying  out  the  operations,  the  above  equation  becomes. 


iexn  00  0 

[Ogxii/21  00  0 

1 0X2i  00  0 

(0x22  0 0 0 

(0gx22/32  0 0 0 

[0gx32/12  0 0 0 

(0x12  0 0 0 

= [0  0 0 0 0 0]T 


Oyll 

Sgyll/21 

Sy21  + 0X21  ( Z21  - 111) 

Oy22  + 0x22  (Z22~  Zll) 
Sgy22/32  + 0gX22/32  (Z22  ~ Zll) 
Sgy32/12  + 0gx32/12  (Z32  - Zll) 
0yl2  + 0x12  (Z12-  Zll) 


Saif  t + 

Ogzii/21  / - 

Sz2i  + 0x2i  (y>2i  - yn)  f + 

0z22 + 0x22  (y>22  - yn)  jT  + 
Sgz22/32  + 0gx22/32  (y22  ~yil)J  + 

8gz32/12  + 0gx32/12  (y32  ~ yil)J  ~ 

Szi2 + Oxn  (y  12-  yii)(T 


Equating  each  component  to  zero  gives  following  three  equations: 

0x11  + 0gxll/21  - 0x21  + 0x22  + 0gx22/32  + 0gx32/12  - 0x12  = 0, 

Syll  + Sgyll/21  ~ Sy21  - 0x21  (Z21  - Zll)  + Sy22  + 0X22  (Z22-  Zll)  + Sgy22/32  + 

+ 0gx22/32(Z22~  Zll)  +Sgy32/12+  0gx32/12(Z32- Zll)  - 5yU  ~ 0X12  (Zl2~  Zll)  = 0, 

8zn  + Sgzll/21  - 0z21  - 0x21  (y21-yil)  + 8z22  + 0x22(y22-  y 11)  + Sgz22/32  + 
0gx22/32  (y22  - yii)  + 0gZ32/i2  + 0gx32/i2  (y32  - yn)  - Szi2  - 0xu  (yn  - yn)  = 0. 


Similar  procedures  would  yield  the  rest  of  the  equations  for  the  remaining  four  loops. 

The  gap  parameters  have  to  be  generated  / constrained  based  on  the  physical  nature  of  mating  between 
the  pair  of  surfaces  and/or  design  functional  requirements.  Converting  a functional  requirement  into 
some  constraint  binding  the  gap  parameters  is  a very  complex  task  and  currently  not  many  generalized 
treatments  have  been  published  (as  far  as  knowledge  of  this  author  goes).  However,  many  people  are 
carrying  out  work  for  this  goal.  In  our  present  work,  we  have  assumed  that  the  mappings  from 
functional  requirements  to  constraints  binding  the  gap  parameters  have  been  carried  out,  and  the 
corresponding  constraint  equations  have  been  obtained. 
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Appendix  - B:  Class  diagrams  of  the  basic  entities  in  the 

synthesis  scheme 
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Appendix  - C:  Artifact  Database  Tables  and  Relations 


TableART 

Purpose  Main  Artifact  Master  File.  Serves  as  the  Artifact  Library. 

Both  compound  and  Primitives  are  stored  here. 
lndexART_NAM  + ART_CAT 
Link- 

StructureVariable  Name  Type  Len  Dec  Details/  Purpose/  Remarks 


ART  NAM 

Char 

10 

QArtifact  Name 

ART  CAT 

Char 

5 

OCategory 

ART  TYP 

Char 

5 

OType 

PRIMITIVE 

Logical 

1 

0 Is  it  a primitive? 

DWGREF 

Char 

15 

0 Drawing  Reference 

LCS  NODE 

Num 

3 

0 Local  Coordinate  system 

CENTER 

Char 

3 

OCenter  Node  (if  any) 

TableARTLOOP 

Stores  Loops  and 

Paths  that  are 

generated  for  each  artifact  for  the  connectivity 

Purposegraph 

IndexART  NAM  + LOOP  +LOOP 

EL 

LinkART  -->  ARTLOOP  (one  to  many) 

StructureVariable  Name 

Type 

Len  Dec  Details/  Purpose/  Remarks 

ART  NAM 

Char 

10 

0 Artifact  name 

LOOP 

Char 

3 

OLoop  counter 

LOOP  EL 

Num 

3 

OLoop  sub-counter 

LOOPS 

Char 

254 

OLoop  String 
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TableARTTREE 

PurposeThe  artifact  link  definition. 

IndexART  NAM  + LOOP  + FROM  ART  + FROM  FACE 


LinkART  -->  ARTTREE  (one  to  many) 

Structure 

Variable  Name  Type  Len  Dec  Details/  Purpose/  Remarks 

ART  NAM 

Char 

10 

OArtifact  name 

LOOP 

Char 

3 

OLoop  number 

FROM  ART 

Char 

10 

OFrom  part 

FROM  FACE 

Char 

3 

OFrom  feature 

TO  ART 

Char 

10 

OTo  part 

TO  FACE 

Char 

3 

OTo  feature 

LOOP  EL 

Char 

3 

OLoop  element 

Table  ARTEQN 

PurposeGenerated  list  of  all  variables  as  used  in  the  constraints  equations  are  kept  here 
This  table  serves  as  a link  between  the  variable  names  and  its  physical  position 
lndexART_NAM+VAR_TYP+FROM_ART+FROM_FACE+TO_ART+TO_FACE+DOF 
LinkART  ->  ARTEQN  (one  to  many) 

; Details/  Purpose/  Remarks 
OArtifact  name 
OVariable  type 
OFrom  part 
OFrom  feature 
OTo  Part 
OTo  feature 
ODegrees  of  freedom 
OVariable  name 
OSerial  number 

0 Equivalent  array  element  name:  z(sln) 
5Value  (computed/specified) 

0 Known  Value? 


TableARTF2N 


Structure  Variable  Name  Type  Len  De 

ART_NAM  Char  10 

VAR_TYP  Char  1 

FROM_ART  Char  10 

FROM_FACE  Char  3 

TO_ART  Char  10 

TO_FACE  Char  3 

DOF  Char  1 

VAR_NAM  Char  15 

VAR_SLN  Num  4 

VAR_XNAM  Char  7 

VAR_VALU  Num  10 

KNOWN_VALU  Logical  1 
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PurposeThis  table  lists  the  feature  to  node  links.  Each  feature  of  each  part  is  defined  as  a 
set  of  connected  nodes. 
lndexART_NAM  + FACE_NO  + NODE_NUM 
LinkART  -->  ARTF2N  (one  to  many) 


Variable  Name 

Type 

Len  Dec  Details/  Purpose/  Rer 

ART  NAM 

Char 

10 

0 Artifact  name 

FACE  NO 

Char 

3 

0 Face/feature  number 

NODE  NUM 

Char 

3 

Onode  number 

DX 

Num 

10 

5delta-x 

DY 

Num 

10 

5delta-y 

DZ 

Num 

10 

5delta-z 

RX 

Num 

10 

5rot-x 

RY 

Num 

10 

5rot-y 

RZ 

Num 

10 

5rot-z 

TableARTFACE 

PurposeTable  defining  each  feature  on  a part 
IndexART  NAM  + FACE  NO 


LinkART  -->  ARTPART  (one  to  many) 

StructureVariable  Name  Type  Len  Dec  Details/  Purpose/  Remarks 

ART  NAM 

C 

10 

OArtifact  name 

FACE  NO 

C 

3 

0 Face/Feature  number 

FACE  TYPE 

C 

3 

0 Feature  type 

CENTER 

C 

3 

OCenter  node(if  any) 

DATUM 

C 

5 

ODatum  reference,  if  any 

DOF 

C 

6 

Odegrees  of  freedom 

DX 

N 

10 

5delta-x 

DY 

N 

10 

5delta-y 

DZ 

N 

10 

5delta-z 

RX 

N 

10 

5rot-x 

RY 

N 

10 

5rot-y 

RZ 

N 

10 

5rot-z 

DCX 

N 

8 

3orientation  direction  cosine-x 

DCY 

N 

8 

3orientation  direction  cosine-y 

DCZ 

N 

8 

3 orientation  direction  cosine-z 

TableARTPART 

PurposeTable  of  parts  of  each  compound  artifact 
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lndexART_NAM  + ART_REQD 
LinkART  -->  ARTPART  (one  to  many) 


StructureVariable  Name 

Type 

Len  1 

Dec  Details/  Purpose/  Remarks 

ART  NAM 

C 

10 

OArtifact  name 

ART  REQD 

C 

10 

0 Requires  (artifact  reference) 

TERMINAL 

L 

1 

01s  this  a terminal? 

CENTER 

C 

3 

OCenter  (node),  if  any 

DOF 

C 

6 

ODegrees  of  freedom  six  'T  or 

DX 

N 

10 

5delta-x 

DY 

N 

10 

5delta-y 

DZ 

N 

10 

5delta-z 

RX 

N 

10 

5rot-x 

RY 

N 

10 

5rot-y 

RZ 

N 

10 

5rot-z 

DCX 

N 

8 

5orientation  direction  cosine-x 

DCY 

N 

8 

5 orientation  direction  cosine-y 

DCZ 

N 

8 

5orientation  direction  cosine-z 

Table  ARTNODE 


PurposeNodal  coordinates  of  all  nodes  on  a part 
lndexART_NAM  + NODE 
LinkART  -->  ARTNODE  (one  to  many) 


StructureVariable  Name 

Type 

Len 

Dec  Details/  Purpose/  Remarks 

ART  NAM 

C 

10 

OArtifact  name 

NODE 

C 

3 

ONode  Number 

CX 

N 

10 

3x-coord 

CY 

N 

10 

3y-coord 

CZ 

N 

10 

3z-coord 

TableARTIO 

PurposeTable  of  I/O  (Input/Output)  parameters  for  a part  / artifact 
IndexART  NAM  + IONAME 
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LinkART  -->  ARTNODE  (one  to  many) 


Structure  Variable  Name 

Type 

Len  Dec  Details/  Purpose/  Remarks 

ART  NAM 

C 

10 

QArtifact  name 

IOTYPE 

C 

5 

01/0  Type 

IONAME 

C 

10 

01/0  Name 

IOC  AT 

C 

5 

0 Category 

iOUNIT 

C 

5 

OUnit 

IOVALU 

N 

10 

3Value 
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Appendix  - D:  Tolerance  Synthesis  Example  Model 


In  this  section  we  present  preliminary  results  for  an  example  artifact  with  three  blocks  (consisting  of 
parts  1,  2,  3)  as  shown  in  Figure  29.  This  artifact  as  represented  in  the  program  is  shown  along  with 
local  coordinate  systems  for  each  feature  in  Figure  30. 

The  program  first  generates  the  artifact  tree  from  the  input  configuration  (Figure  31)  and  then  all 
independent  paths/loops  are  generated  in  an  exhaustive  manner  (Figure  32).  Each  path/loop  is  used  to 
generate  assemblability  constraints.  As  has  been  shown  earlier  these  assemblability  constraints  are  in 
the  form  KX  = b.  For  ease  of  handling  these  constraints  in  Matlab©,  a sparse  matrix  representation  has 
been  used  for  the  ^-matrix.  This  implementation  has  a built-in  procedure  to  generate  cost  functions 
based  on  deviation  parameters.  The  cost  function  generated  for  this  example  is  shown  in  Figure  33. 

For  solving  the  optimization  (minimization)  problem  as  formulated  in  this  program,  we  are  using 
fmincon,  a nonlinear  constrained  optimization  facility  in  Matlab©.  The  program  generates  a matlab 
file  for  running  the  optimization  in  Matlab©.  Results  of  the  optimization  run  are  depicted  in  Figure  34 
as  a distorted  shape  of  the  original  (nominal)  artifact  shown  in  Figure  30. 

These  deviation  parameters  signify  the  amount  permissible  variations  of  each  critical  feature  of  the 
assembly  within  the  prescribed  ranges  of  assembly  and  functional  constraints.  In  order  to  arrive  at  a 
specific  tolerance  specification  we  need  to  carry  out  following  steps: 

o Identify  candidate  features  for  tolerancing, 
o Determine  datum  reference  frames  (DRF), 
o Decide  tolerance  types, 

o Generate  Tolerance  Specification  from  Deviation  to  Tolerance  Mappings  for  the  corresponding  features. 

This  work  is  in  progress  and  we  will  report  actual  mapping  to  a particular  tolerance  specification  in  our 
next  NIST  publication. 


Figure  29:  Block  Artifact 
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Figure  30:  Nominal  shape  with  local  coordinate  systems  (LCS) 
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Figure  31:  Feature  to  feature  connectivity  diagram 
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Figure  32:  Artifact  Control  main  Screen 
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Figure  33:  Sample  cost  of  Manufacturing  Function 
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Figure  34:  Optimized  configuration 
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